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SYNOPSIS 


The various methods of vibration control aim at modifying i 
source or the system or the transmission path from the source 
the system. In vibration isolation, a resilient element 
inserted in the path of vibration transmission from the source 
the system. The simplest idealisation of vibration isolat. 
system is a single degree-of-freedom model in which a mass 
constrained to move only in one direction. The isolator is assu; 
to be massless with linear stiffness and linear damping, 
foundation is considered as ideally rigid. The performa 
characteristics of isolators under the assumption of linear 
have been widely reported in the literature. 

Experiments conducted on real-life isolators such 
wire-rope vibration isolation systems used in space applicatic 
rubber and pneumatic suspensions used in automobile and ot 
applications indicate the need to include nonlinearity in both 
stiffness and damping. In the present thesis, damping foi 
proportional to pth-power of velocity (p=l, 2, and 3, represenl 
viscous, orifice and cubic damping) and combined Coulomb (p=0) 
viscous type are considered. Hard/soft as well 
symraetric/asymmetric restoring forces are included. A brief re^ 
of previous works of Den Hartog, Ruzicka and other researchers 
isolators with linear restoring force and nonlinear damping mo< 
is provided in the first chapter of the thesis. 

The primary objective of this thesis is to analyse the 
of nonlinearity in both stiffness and damping on the perform 



of single degree-of-freedom, passive, vibration isolation systems. 
The governing equation of the isolation system is given by a 
generalised buffing's equation with harmonic excitation. Both 
force and base excitation problems are considered. The method of 
harmonic balance is used to obtain the transmissibility 
expressions. These results are presented in chapter 2. The linear 
stability analysis of the harmonic solution can be carried out 
using the usual variational technique for only odd values of p. 
One encounters analytical difficulties for other values of p. A 
heuristic criterion (based on the notion of equivalent linear 
damping) is proposed to extend the stability analysis which works 
for all values of p (except p=0, i.e., Coulomb damping). Even 

higher order unstable zones can be obtained by this method. 

Today, the role of damping in nonlinear vibration poses a new 
challenge to the theoretician. Even the notions like an attractor, 
basins of attraction etc., are due to damping in the system. While 
tremendous progress has been achieved in understanding various 
features of nonlinear vibrating systems (such as bifurcation 
structure and chaos) , the role of damping and strictly dissipative 
nonlinear damping in particular (as opposed to self-excited 
oscillators like the van der Pol system) has not been addressed to 
in the literature. Therefore, another objective of the present 
thesis is to understand the role of nonlinear damping on the 
bifurcations and chaotic phenomena occurring in harmonically 
excited Buffing's oscillator with pth-power velocity dependent 
damping. In chapter 3, numerical simulation results are presented 
to show that the bifurcation structure and routes to chaos 
inS^nSltlVG to tlio choico of tVlP. ■nonl Of 


are 


In chapter 4, isolation systems with combined Coulomb ^ 

viscous damping and cubic, hard nonlinear restoring force e 

analysed. The method of harmonic balance is used again to obte 

the steady-state response and transraissibility indices with 1 

assumption that the motion is continuous without any stop. 

anomalous jump in the response, similar to the one obtained 

earlier studies on soft systems, is observed when the isolat 
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system is subjected to a base excitation. The phenomenon 
break-loose frequency is examined in detail to bypass the j 
phenomena in the response. 

Numerical simulations on the system discussed in chaptei 
are carried out in chapter 5 to analyse the role of Coul 
damping in bifurcation structure and chaos. It is shown that fo 
base-excited system, the inclusion of a Coulomb damper witl 
suitable break-loose frequency can suppress the second 
resonances and chaotic motion. However, for a force-exci 
system, the introduction of Coulomb damping does not alter 
bifurcation structure. 

Soft buffing's oscillators (possessing a saddle point in 
unforced situation) such as, single- and double-well potent 
oscillators represent a variety of engineering and phys' 
problems. These systems are widely studied in literature as s 
analytical work can be carried out due to the existence of 
homoclinic orbits. The role of nonlinear dissipation on 
bifurcation structure in these systems is studied in chaptei 
The Melinkov criterion and an analytical criterion for 
period-doubling are obtained in the presence of combined vis 
and cubic damping. 


xii 

Some preliminary experimental results are presented in 
chapter 7 for a soft isolator (made in the form of a hollow rubber 
tube) . The transmissibility characteristics of this system 
obtained experimentally indicate that the soft characteristics of 
the hollow rubber tube are beneficial for effective vibration 
isolation. 

The main conclusions of the thesis are: 

(i) A nonlinear isolator with a soft characteristic is superior 
to that with a hard characteristic. 

(ii) For a base-excited system, a suitable choice of the Coulomb 
damping coefficient enables one to maintain the transmissibility 
at or less than unity over the entire frequency range. 

(iii) Both the period-doubling and intermittency routes to chaos 
occur for all values of the damping exponent p, thus preserving 
the bifurcation structure of a viscously damped, Duffing 's.i 
oscillator. It is shown that strictly dissipative nonlinear 
damping can be effectively used for controlling chaos. 

(iv) It has been demonstrated by numerical simulation that in a 
base-excited system, the instabilities associated with the 
secondary resonances and chaotic motion, encountered in the low 

frequency regime, can be eliminated by using a suitably designed 
friction damper. 

(V) For a double-well potential oscillator, with the same value 
of the damping coefficient, the threshold values of forcing 
required to cause period-doubling bifurcation are seen to decrease 
with increasing value of the damping exponent. This feature is 

Just opposite to what has been observed in the case of a hard, 
Duffing's oscillator. 



(vi) Experimental results indicate that a suitable combination 
a soft nonlinear spring and the optimum value of the damping 
yield a better design than the linear isolator. 
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CHAPTER I 


INTRODUCTION 


1.1 Introduction 

The quest for safe, reliable and quieter technologic 
increasingly becoming important in the present day v 
Vibration and noise control engineering is an important comj 
of such an endeavour. Light weight vehicles and machines ri. 
at high speeds, performance sensitivity of electronic and o£ 
instruments, design of robust launch vehicles and 

structures, protection of buildings from earthquakes and 
natural calamities give a fresh impetus for rapid developme 
new techniques of vibration control. The awareness of ergt 
aspects in designing shop floors and buildings etc. , to ere 
safe and pleasant working environment devoid of nois< 
vibration is another source of motivation for the progre 
vibration control. 

In many engineering situations, the source of vibrati 
either the inertia forces of the moving parts or the oscil 
foundation/base motion. In such cases, it may not be feasil 
reduce the level of excitation below the acceptable lin 
modify the design so as to bring down the response belc 
desired level. Therefore, it is apparent that secondary meth 
control are necessary to tackle many vibration problems. Vib 
isolation is one such important and widely used method. 

A vibration isolator in its most elementary form i 
considered as a resilient^ member inserted between the sour 
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receiver of vibration. Many different materials arc used in the 
form of isolator. Isolators can be rubber springs, air spring.s and 
coiled or leaf metal springs. Each one of thor.e has some 
advantages and disadvantages. The damping capacity of an isolator 
may originate from various sources. For a bonded-rubber spring, 
energy dissipation occurs due to the inherent laige internal 
friction within the material. With a coiled metal spring, normally 
a viscous dashpot is used in parallel to provide the damping. With 
leaf-springs, interface slip damping is a major contributor 
towards the overall damping. In the case of air-s-pr i ngr, , orifices 
placed in the flow path cause the desired energy disr.ipat ion . Most 
of these damping mechanisms, like viscous. Coulomb (si Lp) or 
orifice, can be adequately treated as special cases of n damping 
force proportional to pth-power of velocity. This, velocity 
dependent model, however, fails to account for the hysterotic type 
of damping force. 

The function of a vibration isolator is primarily to reduce 
the magnitude of the motion transmitted from a vibrating 
foundation to the equipment (referred to as the bas.c excitation 
problem) , or to reduce the magnitude of the force trans.mitted f rom 
the equipment to the foundation (referred to as the force 
excitation problem) . Dependinq on the method of producing the 
controlling force, vibration isolation can be classified as 
passive, active and semiactive. In active vibration isolation, the 
isolator force is servo-controlled that enhances the system 
performance. But these active control schemes are expensive and 
are usually limited to light weight machines, if the isolator 
force is not controlled using a feed-back strategy, then the 



method is referred to as passive vibration isolation. A 
combination of the above two schemes is referred to as seniactive 
vibration isolation. 

The simplest idealisation of a vibration isolation system is 
a single degree-of—freedom model in which a mass is constrained to 
move only in one direction. The isolator is assumed to be uassless 
with linear stiffness and linear damping. The foundation is 
considered as ideally rigid. More complex vibration isolation 
systems can be visualised by relaxing these assumptions, such as 
when (i) the foundation is not ideally rigid, (ii) the isolator 
inertia cannot be ignored, (iii) several masses with intermediate 
elastic elements are involved in the model or (iv) the mass can 
move in any direction producing e.g., heave, roll, yaw and pitch 
modes. With the advent of modal analysis still more complex models 
treating the machine as a structure (no longer a rigid mass) 
connected to the foundation structure by isolators at several 
points are also not uncommon in the literature. But with 
increasing complexity, the problem becomes more specific and it is 
difficult to obtain generalised conclusions. 

Further complications of the vibration isolation system may 
arise out of the type of excitation such as random or impulse 
excitation. The performance characteristics of isolators under the 
assumption of linearity have been widely reported in the 
literature. For such linear systems under harmonic excitations, 
various transfer functions, e.g., displacement or force 
transmissibility are used as the indices of effectiveness of an 
isolator. The various indices are related to one another through 
simple expressions and the choice of a particular index depends on 
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the type of the response variable and applications. Though the 
research on vibration isolation is vast, very little attention is 
paid to provide guidelines to design or select the isolators on a 
rational basis. The general practice of designing of the isolators 
is still on the assumption that all will be well if the natural 
frequency of the isolated system is appreciably lower than the 
frequency of excitation. 

Consideration of the nonlinearities present in most of the 


practical isolators becomes important specially during the design 

and analysis of suspensions of high-speed vehicles, mounts for 

optical and other sensitive instruments to be used in space 

vehicles. Experiments conducted on real-life isolators such as 

''^it'e-rope vibration isolation systems used in space applications, 

rubber and pneumatic suspensions used in automobile and other 

applications indicate the need to include nonlinearity in both the 

stiffness and damping. While, there has been a significant 

progress in understanding the role of nonlinearity in mechanical 

systems from the point of view of geometric theory of nonlinear 

dynamics, these results have not percolated down towards 

designing/ improving the existing isolation systems. In the 

nonlinear theory of vibration isolation, one encounters with the 

problem of defining a suitable performance index for the isolator. 

This is because a harmonic response with^same frequency as that of 

the excitation is not guaranteed. The response may contain 

subharmonics , superharmonics and somet-iTnoc 

ana sometimes the response may even 

be nonperiodic (chaotic). Furthermore, unlike in f, t 

/ uniiKe in a linear system, 

the indices to be used for various types of PVf>ii.t.+.- 

yp ot excitation cannot be 

related through simple expressions and 


consequently force and base 



excitation problems need to be analysed separately. 

The primary objective of this thesis is to analyse the role 
of nonlinearity in both the stiffness and damping on the 
performance of vibration isolation systems. In the present work, 
the isolation system is modelled as a generalised Buffing's 
oscillator with pth-power damping subjected to a harmonic 
excitation. Though the results presented are in the context of 
vibration isolation, emphasis is also laid on understanding the 
role of strictly dissipative nonlinear damping in the general 
framework of the geometric theory of dynamical systems. 

1.2 Review of Previous Work 

The amount of literature available on the general problem of 
vibration isolation is too vast to be considered in detail. In 
this section, the literature pertinent to the present work is 
discussed. A summary of the investigations reporting various 
features of the Buffing's oscillator is also presented. An 
overview of the research work in different areas, connected with 
the present problem, is categorized under the following headings: 

(i) Vibration isolation, 

(ii) Buffing's equation, 

(iii) Systems with a Coulomb damper. 

1.2.1 Vibration isolation 

The performance characteristics of isolators and other 
results obtained prior to 1950 have been summarised in Crede's 
text book [1]. Betails of a number of commercial isolators 
including their hardware constructional features have been 



presented in reference [2]. A recent account of the research 
carried out can be obtained from reference [3]. Researchers have 
used different indices to express the effectiveness of an 
isolation system. Some of these are: (1) Force/displacement 
transmissibility (absolute), (2) Relative transmissibility , (3) 

Response ratio, (4) Isolation effectiveness, (5) Insertion loss, 
and (6) Vibration power flow. Log-Log plots of the above mentioned 
quantities vs. frequency ratio can be found extensively in the 
literature for various values of the system parameters. 

Vibration isolation problems are generally analysed with the 
assumption that the foundation possesses infinite mechanical 
impedance. This assumption, although valid in many practical 
situations, is not applicable when the vibratory system is mounted 
on, say, ship, aircraft or automobile structures. In such a 
situation, the finite impedance of the foundation drastically 
affects the vibration isolation characteristics. Snowdon [4] has 
studied this problem when the isolator is made of rubber like 
materials for rigid as well as non-rigid foundations. The 
frequency dependence of the dynamic modulus and damping factor of 
rubber like materials has been taken into consideration while 
calculating the above mentioned indices. The concept of parallel 
combination of low and high damping rubber for improved 
performance of the isolation system has been proposed. Snowdon has 
used ^response ratio^ to quantify the effect of vibration 
isolation for nonrigid foundations. The nonrigid foundations are 
modelled as simply-supported beams and plates. This analysis has 
been extended to the case when the feet of the machine supported 


on the isolator are nonrigid [ 5 ]. m this reference, 


Snowdon has 
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presented an excellent review of the use of elastomeric materials 
in vibration isolation. Wave effects may be observed at high 
frequencies when the rubber mount dimensions become comparable 
with multiples of the half-wave lengths of the elastic waves 
travelling through the mounting. Alternatively, wave effects may 
be thought of as occurring when the elasticity and the distributed 
mass of the rubber mounting interact at high frequencies. A guide 
to the character of the wave effects in antivibration mountings 
has been obtained by assuming that the 'mountings' obey the simple 
wave equation for the longitudinal vibration of a rod of uniform 
cross section [4]. Both 'long rod' and 'Love' theories have been 
used to analyse the wave effects [5]. Also it has been shown that 
the performance of an isolator is overestimated by ignoring the 
inertia and the resulting wave effects. 

Macinante [6] has considered a two mass model to take care of 
the flexibility of floors and carried out a detailed parametric 
study. In reference [7], the criteria of effective isolation of 
machinery have been derived and the design of isolators complying 
with these criteria has been described. The concept of vibration 
power flow has been explained by Goyder and White [8] . Expressions 
have been derived for vibration power flow to the foundation for 
single point excitation, taking the foundation mobility into 
consideration. Design considerations for isolators and 
foundations, to minimise the transmission of vibration power in 
the structure, are examined. Pinnington and White [9] have also 
investigated the parameters controlling the power transmission 
from machine to the seating structure via vibration isolators. 
Some vibration isolation problems encountered on ships have been 
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considered in reference [10]. 

If added mass can be tolerated, a two stage or compound 
mounting can provide low values of transmissibil ity especially at 
high frequencies. The design of a practical mount based on this 
idea has been reported in references [4, 5] where a complete 
analysis of the compound stage mounting with rigid and nonrigid 
foundations has ben given. High frequency vibration isolation of 
two-stage isolators has been discussed in reference [11]. Some 
general design criteria for non— dissipative vibration isolation 
systems have been given in reference [12] . An analytical procedure 
for the evaluation of the transmissibility of N degreos-of-f reedom 
viscoelastic mountings has been developed [13]/ where the 
transmissibility is expressed through a class of polynomials. The 
use of four-pole parameters which enables one to take the isolator 
inertia into account has been demonstrated in references [5, 14], 
In reference [15] the wave effects in an isolator, consisting of N 
equal masses and N intermediate distributed elastic mounts, was 
analysed by using the four-pole parameters. The resulting 
transmissibility expressions are lengthy polynomials, especially 
for increasing values of N. The coefficients of these polynomials 
have to be evaluated separately for each value of N. In reference 
[16] the notion of a propagation constant has been used profitably 
to render the computational effort independent of the number of 
periodic elements present in the system. The analysis includes 
isolator inertia, frequency dependence of the isolator stiffness 
and damping, and the finite foundation impedance. Periodic 
isolators (antivibration mounts consisting of rubber blocks with 
intermediate steel spacers) are shown to 


possess better high 
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frequency isolation characteristics as compared to a simple one 
with the same overall dimensions. 

Pinnington [17] has derived the expressions for vibration 
power flow when a motor is supported by four isolators on a 
flexible seating system. A detailed dynamic analysis of a heavily 
damped foundation structure designed for the multidirectional 
isolation of a diesel engine has been presented in reference [18]. 
The problem of defining suitable performance indices of an 
isolation system subjected to random excitations has been 
addressed to in reference [19]-. A detailed review of the 
literature on vibration isolation has been compiled in reference 
[20] under the following headings: (i) static and dynamic 
properties of isolator materials, (ii) stress analysis of rubber 
mounts, (iii) modelling of finite impedance foundations, (iv) 
analysis of a single point excitation model, (v) analysis of an 
unidirectional multidegrees-of-freedom model, (vi) analysis of two 
or multipoint excitation model and, (vii) wave effects in 
isolators. 

A better model of a number of practical isolators should 
include the nonlinear characteristics inherent in the spring and 
the damper. A study of the dynamic characteristics of a wire rope 
isolation system (constructed with helical strands) used in 
numerous space and military applications has been presented in 
reference [21]. A semi-empirical model having non-linear 
stiffness, power-law damping and variable Coulomb friction damping 
has been developed, and the results have been compared with the 
experimental data. Ruzicka and Derby [22] have summarised the 
results of scores of individual investigations on the problem of 
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vibration isolation of a mass with single degree-of -freedom. 
Nonlinear damping models considered in this work include linear 
dampers with clearance, Coulomb friction dampers, hysteretic 
dampers, and dampers whose damping force is proportional to a 
general power of velocity across the damper. Several indices of 
the performance of the isolation system with such nonlinear 
damping models and linear stiffness characteristics have been 
presented in a form which is immensely useful to a design 
engineer. 

Kirk [23] has presented a study of the influence of nonlinear 
spring stiffness-characteristics on the effctiveness of vibration 
isolators with linear damping, subjected to stationary, random 
white noise ground acceleration. The results obtained with the 
cubic hard spring, the cubic soft spring and the tangent spring 
have been compiled as design charts. Isolators having symmetric 
and asymmetric nonlinear restoring force characteristics with 
linear damping have been investigated in reference [24], where the 
'jump' phenomenon observed in the response is explained through 
the bifurcation set and catastrophe theory. The results are 
compared with those obtained from an experiment conducted with a 
pneumatic spring of constant area diaphragm. The application of 
catastrophe theory, rather than the study of isolation 
characteristics, seems to be the main objective of the work 
presented in reference [24]. 

A two degrees-of-freedom model incorporating nonlinear 
stiffness and damping characteristics is reported by Metwali [25]. 
Complex motion of the isolated mass at high excitations, taking 
into account the effects of the nonlinear material behaviour on 
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transmissibility , has been studied and bifurcations have been 
observed [26]. Beatty has considered the finite amplitude, 
periodic motion of a body supported on rubber shear springs [27, 
28]. This work has been extended in reference [29]. The effect of 
the geometric and material nonlinearities of the elastomeric 
mounts on the transmissibility of a simple system has been studied 
experimentally by Harris [30, 31]. The types of nonlinearities 
encountered in mounts of various geometries and modes of loading 
have been obtained. It has been pointed out that these nonlinear 
aspects of rubber mounts have been relatively unexplored for 
antivibration devices. Also the role of nonlinearity in designing 
better suspensions for automobiles, nonwheeled ground vehicles and 
vibration control problems in space applications has been 
discussed. The standard approach to vibration isolation analysis 
represents the source, the isolator and the receiver by their 
mobilities. The choice of the appropriate description for each 
element is not arbitrary, specifically when the isolation system 
is nonlinear or when the system has noise. At least one of the 
system components must be represented by its impedance [32]. An 
interesting discussion of the challenges encountered in the low 
frequency vibration isolation, required to carry out microgravity 
science experiments has been presented in reference [33]. In 
reference [34], design and control aspects of microgravity 
isolation mounts have been detailed. The steady-state, harmonic 
response and the transmissibility indices of nonlinear isolators 
have been discussed in reference [35]. 
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1.2.2 Duffing’ s equation 

In the present work, the governing equation of motion of the 
isolation system is given by a generalised Duffing’s oscillator 
with nonlinear damping. Hence, a summary of the recent 
investigations, carried out on various features exhibited by 
buffing's oscillator and the role of damping, is presented in this 
section, buffing's equation is one of the paradigm cases of 
nonlinear vibrating systems. Classical results on buffing's 
equation can be obtained from the standard text books [36-42]. A 
modern treatment of buffing's oscillator from the geometrical 
theory of dynamical systems has been summarised in references 
[43-50] . 

buffing's oscillators are employed as models of various 
physical and engineering situations such as Josephson junctions, 
optical bistability, plasma oscillations, buckled beam, ship 
dynamics, vibration isolators and electrical circuits, etc., [35, 
43-52]. Since Ueda's [53-55] work on buffing's equation, it is 
well-known that nonlinear systems under consideration can exhibit 
chaotic responses (under harmonic excitation) in certain parameter 
regimes. Such chaotic responses are studied via phase plots, 
Poincare maps, power spectrum, Lyapunov exponents etc. [43-45]; all 
obtained by using numerical methods. Single degree-of-freedom 
systems with different types of nonlinear restoring force and a 
linear damping force have been studied by various workers. With a 
symmetric nonlinear restoring force, the existence of a 
symmetry-breaking precursor and subsequent period-doubling route 
to chaos have been explained through the stability analysis using 
Hill's equation [56]. 
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Crutchfield and Huberman [57, 58] considered the role of 
fluctuations on the onset of the period-doubling route to chaos. 
They showed that the structure of the strange attractor is very 
stable even under the influence of large fluctuating forces and 
that the role of noise is to introduce a symmetric gap in the 
deterministic bifurcation sequence. The effects of noisy 
precursors and the so-called virtual Hopf bifurcation on the 
period-doubling have been studied by Wiesenfeld [59-61]. 

Parlitz and Lauterborn [62] examined the connection between 
the bifurcation set and resonance structure of a forced Duffing's 
oscillator. In a series of papers [63-66] they employed the notion 
of 'Torsion Number' to classify resonances and discussed the 
universal superstructure in driven, dissipative nonlinear 
oscillators. These results indicate various scenarios in the 
bifurcation structure of Duffing's equation viz., hysteresis, 
co-existence of attractors, symmetry-breaking, existence of dual 
responses, period-doubling route to chaos. The work has also been 
extended to construct two-dimensional maps capturing similar 
features [66]. Fang and Dowell [67] also observed dual responses 
in Buffing's equation. The existence of intermittency transition 
to chaos in a symmetric potential well is explained in reference 
[68]. Floquet theory and harmonic balance method have been used to 
calculate the transition boundaries. In a series of papers, 
Szemplinska-Stupnicka [69-71] examined the role of classical 
methods in determining the subharmonic and chaotic zones in the 
amplitude-frequency plane. Also it has been demonstrated that the 
classical stability analysis using Mathieu and Hill's equations 
can capture the period-doubling observed in numerical simulations. 
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Several authors [72-80] obtained an extra jump in the 
harmonic response of soft Duffing oscillators indicating a break 
in the resonance curve. Miles [75, 76] referred to this jump as an 
'anomalous jump^ and conjectured that this break in the resonance 
curve for the oscillator is a necessary antecedent to 
symmetry-breaking . 

In all the literature cited so far, the phenomenological 
model of the dissipative force has been assumed to be linear. In 
this context it is appropriate to recall what Pippard [81] has 
said “There is something of a tendency among physicists to try to 
reduce everything to linearity...... reality may not always conform 

to what we might wish, rather more so with the damping forces than 
with the restoring force in small-amplitude vibrations.". In the 
context of his discussion of the experimental results obtained by 
Wraight [82], Pippard also suggested that the nonlinear damping 
curves may give a measure of the pinning and frictional forces 
involved when the magnetic fields penetrate and move within a 
superconductor. The role of dissipation in Josephson junctions has 
received special attention in recent years both from the 
experimental [83] and theoretical [84] points of view. Landauer 
[85] presented a nice discussion of ensuing philosophical debate 
on the origin of dissipation, fluctuations and irreversibility. 
The modern geometrical theory of the dissipative dynamical systems 
introduced several new notions like 'strange attractors', 'basins 
of attraction' etc, [86]. These concepts, being essentially the 
facets of dissipation in dynamics, seem to inherit all the 
associated problems. Milnor [87] gave a detailed account of the 
conceptual difficulties for giving a precise mathematical 
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definition of an attractor. 

A good account of the role of damping in vibrations from the 
engineering point of view has been given by Crandall [88]. He 

mentioned that damping is of great relevance as it decides the 
border of stability and instability. Consideration of nonlinear 
damping models is necessary in several engineering applications 
such as the "roll damping" in ship dynamics [51, 52], vibration 

Isolators [21, 35] and the drag forces encountered in flow induced 

vibration problems [41, 44]. The importance of nonlinear damping 
in engineering stems from the fact that it can be used as an 

effective passive control strategy to suppress various 
instabilities . With these applications in mind a study of the 
role of nonlinear damping in the dynamics of Buffing's oscillator 
has been carried out in references [89, 90]. 

1.2.3 Systems with Coulomb damping 

Detailed performance characteristics of a vibration isolation 

system with directly as well as elastically coupled Coulomb damper 

have been presented in reference [22]. The notion of equivalent 

viscous damping was employed to obtain expressions for various 

transmissibility indices. The importance of the break-loose 

frequency has been explained and the variations of the high 

» 

frequency attenuation rate and resonance transmissibility with 
various values of damping have been discussed. An exact solution 
for the symmetric, steady-state response of a vibrating system 
with combined Coulomb and viscous friction has been given by Den 
Hartog [91]. He considered the motion with and without stops and 
examined the range of validity of the solution obtained by 
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assuming continuous motion. An analysis of similar vibrating 
systems, with combined Coulomb and viscous friction, subjected to 
base excitation can be found in references [92, 93]. However, in 
all the references mentioned above, the analysis has been 
restricted to motion with atroost two stops per cycle. An extension 
of Den Hartog's work to obtain response with multiple stops per 
cycle is carried out in reference [94]. Shaw [95] extended the 
work of Den Hartog by considering a situation where the static 
coefficient of friction is different from the kinetic coefficient. 
He also considered the effects of negative viscous damping 
mechanisms which may arise out of aerodynamic forces acting on 
turbine blades. The stability analysis was carried out by using 
bifurcation theory and the possibility of existence of aperiodic 
motions has been pointed out. It has been shown that the symmetric 
motion with two stops per period can be unstable and that pairs of 
unsymmetric motions are generated at the bifurcation points- It 
has been reported that beating type motions may occur when the 
equivalent viscous damping coefficient is negative. Some transient 
motions have also been analysed. 

Studies of vibrating systems with various types of friction 
mechanisms are also important in the context of protection of 
structures from earthquakes by sliding isolation systems as 
reported in reference [96]. In this work, the transient and steady 
state responses of a force-excited system in the presence of 
various types of damping forces have been reported. The motion 
under consideration consisted of a number of stops per cycle. It 
was concluded that with increasing viscous damping, the number of 
stops per cycle of motion reduces and may even come down to zero. 
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In a series of papers, Dowell and his co-workers [97-99] have 
investigated vibrating systems with friction damper. The method of 
harmonic balance was used in these works. The existence of a 
strange attractor in a single degree-of-freedom nonlinear 
oscillator with friction is reported by Awrejcewicz [100], This 
work has been extended in reference [101], where a detailed 
analysis to determine the parameter zones of chaotic motion has 
been carried out through numerical evaluation of.Pioncare maps and 
Lyapunov exponents. A two degrees-of- freedom, self-excited 
oscillator with friction has been considered in references [102, 
103], where the stability of the equilibrium points received 
special attention. Different types of stick-slip transition zones, 
numerical investigation of periodic and chaotic motions and some 
results on the transient behaviour have also been presented. 
Chaotic motion in self-excited systems with friction damping has 
also been dealt with in reference [104]. 

In reference [105], a hard, Duffing-type (i.e.,with cubic 
nonlinear restoring force) isolation system with combined Coulomb 
and viscous damping subjected to a harmonic excitation has been 
investigated. Both force and base excitation models are presented. 
The role of subharmonics and chaotic motion in the context of 
vibration isolation using Coulomb dampers has been addressed to in 
reference [106]. 

1.3 Objectives and Scope of the Present Work 

The primary objective of the present thesis is to study the 
performance characteristics of a single degree-of-freedom passive 
vibration isolation system with nonlinearity in both spring and 
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damper. The governing equation of the isolation system is given by 
a generalised Duffing's equation with harmonic excitation. A 
pth-power velocity dependent damping force and hard/soft as well 
as symmetric/asymmetric restoring forces are included. Both force 
and base excitation problems are considered. The method of 
harmonic balance is used to obtain the transmissibility 
expressions and the results are presented in a format which can be 
readily used by the designer of isolation systems. 

In 1971, Crandall [22] mentioned in his foreword to a 
monograph on vibration isolation that vibration theory was 
essentially complete except for a realistic treatment of damping. 
He pointed out that though the theoretical advances helped to gain 
a better understanding of the coupling of various modes due to the 
damping in multidegrees-of-freedom linear systems, the problems of 
damping still remain a principal challenge in vibration 
technology. Today, the role of damping in nonlinear vibration 
poses a new challenge to the theoretician. As already mentioned in 
section 1.2.2., even the notions like an attractor, basins of 
attraction etc., are due to damping in the system. While the 
literature survey indicates that tremendous progress has been 
achieved in understanding various features of nonlinear vibrating 
systems, the role of damping and strictly dissipative nonlinear 
damping in particular (as opposed to self-excited oscillators like 
Van der Pol system) has not been addressed to in the literature. 
Therefore, attention is paid to understand the role of nonlinear 
damping on the bifurcations and chaotic phenomena occurring in 
harmonically excited Buffing's equation. 

Recently, exploitation of nonlinearity to obtain better 
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designs is increasingly becoming important [107]. Several attempts 
are being made to devise control schemes to quench the resulting 
bifurcations and chaos [108, 109]. In the present work, an effort 

is also made to understand the role of strictly dissipative 
nonlinear damping as a passive control strategy to suppress 
various instabilities occurring in nonlinear isolation systems. 
The organisation of the thesis is detailed below. 

In chapter 2, a wide class of single degree-of-freedom 
vibration isolation systems subjected to both force and base 
excitations are considered. The nonlinear restoring and damping 
forces are assumed to be given by power laws. Symmetric/asymmetric 
as well as hard/soft spring models are considered. The soft 
systems are assumed to possess no saddle point in the unforced 
case. The method of harmonic balance is used to obtain the 
transmissibility expressions and a parametric study is reported to 
evaluate the performance characteristics. A heuristic procedure 
(using the notion of equivalent viscous damping coefficient) , for 
analysing the stability of the harmonic solution of a hard, ^ 
Duffing 's oscillator with pth-power velocity dependent damping 
(p20, i.e., excluding Coulomb damping) is presented. 

In chapter 3, a special case of the systems considered in 
chapter 2, i.e., Duffing's oscillator with pth-power velocity 

dependent damper subjected to a harmonic force excitation is 
studied using numerical simulation. Phase plots and Poincare maps 
are obtained to study the role of damping on bifurcations and 
chaos. It is shown that the bifurcation structure and routes to 
chaos are insensitive to the choice of the nonlinear damping 


model . 
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In chapter 4, isolation systems with combined Coulomb and 
viscous damping and cubic, hard nonlinear restoring force are 
analysed. The method of harmonic balance is used again to obtain 
the steady-state response and transmissibility indices with the 


assumption that the motion is continuous without any stop. An 
anomalous jump in the response, similar to the one obtained in 
earlier studies on soft systems, is observed when the isolation 
system is subjected to a base excitation. The phenomenon of 
break-loose frequency is examined in detail to bypass the jump 
phenomena in the response. The obtained results extend the 
previous works of Den Hartog and Ruzicka who considered a linear 
restoring element. Numerical simulations are carried out to 
ascertain the status of anomalous jump and the validity of the 
results obtained by the method of harmonic balance. It is pointed 
out that unlike in the case of a soft system, this anomalous jump 
does not seem to have any bearing on symmetry-breaking and 
subsequent period-doubling. 


In chapter 5, it is shown that for a base-excited system, the 
inclusion of a Coulomb damper with a suitable break-loose 
frequency can suppress the secondary resonances and chaotic 
motion. However, for a force-excited system, the introduction of 
Coulomb damping does not alter the bifurcation structure. 


In chapter 6, the effect of nonlinear dissipation on the 
response and bifurcations of soft Duffing's oscillator (possessing 
a saddle point in the unforced situation) is considered. A 
parametric study is reported to indicate the dependence of 
threshold values of the parameters, at which bifurcations occur, 
on the damping index and the damping coefficient. The Melinkov 
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criterion and an analytical criterion for the period-doubling are 
obtained in the presence of combined viscous and cubic damping. 

In chapter 7, experimental results are presented for a 
nonlinear vibration isolator made in the form of a hollow rubber 
tube. 



CHAPTER 2 


HARMONIC RESPONSES OF NONLINEAR VIBRATION ISOLATORS 
WITH pth-POWER DAMPING 


2.1 Introduction 

The performance characteristics of isolators under the 
assumption of linearity have been widely reported in the 
literature [2, 4]. For such linear systems under harmonic 
excitations, various transfer functions, e.g., displacement or 
force transmissibility are used as the indices of effectiveness of 
an isolator. These various indices are related to one another 
through simple expressions and the choice of a particular index 
depends on the type of response variable and applications. 

In the nonlinear theory of vibration isolation, one 
encounters with the problem of defining a suitable performance 
index for the isolator. This is because a harmonic response with 
the same frequency as that of the excitation is not guaranteed. 
The response may contain subharmonics and superharmonics, and 
sometimes the response may even be non-periodic (chaotic) [ 43 3 - If 
other harmonics are present, a working index of isolation 
effectiveness may be defined as the ratio of the r.m.s. values of 
the response and the excitation. In the case of a chaotic 
response, the ratio of the power spectral densities of the 
response and the excitation may be used. These indices provide an 
overall measure of transmissibility [19]. However, even these 
working indices are devoid of any information insofar as the 
frequency content of the transmitted vibrations is concerned. The 
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information regarding frequency content is sometimes essential, 
especially when the foundation is of finite impedance. 
Furthermore, unlike in a linear system, the indices to be used for 
various response variables cannot be related through simple 
expressions and each case needs to be analysed separately. In this 
chapter, however, the standard index of transmissibility is 
retained since only the harmonic response is considered. A minor 
modification, in case the restoring force is asymmetric, has been 
made to take care of the presence of a frequency independent term 
in the response. The method of harmonic balance is used to obtain 
the steady-state, harmonic response and transmissibility under 
both force and base excitations. 

In the present work, the isolator is modelled as one which 
provides nonlinear restoring and damping forces. The usual pth- 
power model has been used for damping which encompasses several 
special cases like Coulomb, linear, and cpaadratic damping models 
with p = 0, 1, 2, respectively. The case of Coulomb damping (i.e., 
p = 0) shows many special features. Therefore, the results with 
p=0 are not included here and will be reported in a separate 
chapter. However, both symmetric and asymmetric nonlinear 
restoring forces are investigated. These models for isolators are 
in line with those suggested in references [2, 24, 25]. Force and 
base excitations are considered separately in each case. With 
harmonic force excitation on a rigid mass (the source) , the 
foundation (the receiver) has been assumed to be fixed (i.e., of 
infinite impedance) . The jump phenomenon in the response similar 
to that of Puffing's equation is obtained. The effects of system 
parameters on the "jump", resonance transmissibility and 
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high-f rquency attenuation rate are presented. 

2.2 Modelling of Isolation System and Governing Equations of Motion 
As stated earlier, two different types of excitations are 
investigated. In the first one (Figure 2.1), a rigid mass M is 
acted upon by a harmonic force f cos wt, and this mass is isolated 
from a fixed foundation through an isolator. Hereafter, this 
situation will be referred to as "force excitation". The second 
type of excitation, hereafter, referred to as "base excitation" is 
explained in Figure 2.2. Here the foundation is moving with a 
harmonic displacement y = yQ cos wt. The governing equations of 
motion for these two excitations with symmetric and asymmetric 
restoring force characteristics of the isolator are obtained 
separately as detailed below. 

2.2.1 CASE I: Isolation system with symmetric restoring force and 
subjected to force excitation 

The governing equation of motion for the system under 
consideration can be written as [25], 

Mx^'+Cx' jx' |^^ + Kx|x|'^^ = f cos ut (2.1) 

wheje a prime denotes differentiation with respect to time t, w 
is the frequency of excitation, C and K are the constants of 
proportionality for damping and restoring forces, respectively. 
The exponents p and q characterise the nature of nonlinearities 
present in the damping and stiffness, respectively, of the 
isolator. This model can be considered as an idealisation of a 
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variabls disphragm pnsuinatic suspension 3s suggested in reference 
[25]. With q=l, equation (2.1) reduces to that considered by 
Ruzicka [22] and with q=3, we get a generalised Duffing 's equation 
with pth-power damping model. It should be mentioned that a 
special case of equation (2.1) with p=l, q=2 has been studied by 
Hayashi [37]. 

Let us define the following non-dimensional parameters: 


T = Wq t, n = W/Wq, X = X/Xq, 

X = xVXqWq, X = x'VXq C = C xP'^ w^“V( 2 M) with, 

Xq = (f/K)^/'5, Uq = [K x^"^ /M]^/^ 

where a dot at the top denotes differentiation with respect to the 
non-dimensional time x. Using the above parameters, equation (2.1) 
is rendered in the following nondimensional form: 

X + 2^X jX|P^ + X 1x1*^^ = cos Qx (2.i 

2.2.2 CASE II: Isolation system with symmetric restoring force and 
subjected to base excitation 

From Figure 2.2, the relative displacement 5 of the isolator 
is given by, 5 — x y, where x is the response and y is the base 

excitation. The governing equation of motion for this system can 
be written as 


M 6"'+ c S' |5' 


I P-1 


+ K 5|5 


-M Y" 


My^ w'^cos wt. 


(2.3) 
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Defining the following parameters, 

T = w^t, Q = w/Uq, a = S/Yq, 

A = S'/ YqWq, a = S'VYq and C = C yg'^ wg"2/(2 M) , with 
Oq = [K 

and substituting in equation (2.3) we get the following 
nondimensional equation; 

A + 2 C A 1A|P"^ + A cos Qr (2.4) 

2.2.3 CASE III; Isolation system with asymmetric restoring force 
and subjected to force excitation 

Let us consider an asymmetric restoring force given by + 

2 3 ... 

K 2 X + K^x . The governing equation of motion can now be written as 

M x" + C x' I x' I*’ * + K^x + K 2 X^ + K 2 X*= f cos Lot (2.5) 

which represents the motion of a mass on a constant-area diaphragm 
pneumatic isolator as shown in references [2,24]. This equation 
with p=l has also been considered in references [37,69]. 

Defining the following nondimensional parameters; 

T = Wq t, n = CJ/Wq, X = X/Xq, 
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X = xVx.u., X = x'Vx. c = c xP ^ cjP V(2 M), 


0 0 


= K^Xq/(.Mu>1) = K 2 Xq/K^, C 2 = K 3 xp K 3 , with 


Xq = f/K^, Wq = [K^/M]""", 


and substituting these in equation (2.5), one obtains 


X + 2 C X |X|P“^ + X + e^X^+ cos Qx . 


( 2 . 6 ) 


Using a transformation X = X + c^/( 3 c 2 ), equation (2.6) can 
be changed to 


X + 2 C X X 




+ [1 - cp (3 c^)] X + e X’ 


cos Qz + c^/ (3 c^) - 2 cj /(27 c2) . 


(2.7) 


Since the presence of the linear term does not change the 
result qualitatively [37], for the sake of simplicity is chosen 
such that the linear term goes to zero. Thus, upon choosing = 
/"3e^, equation (2.7) reduces to 


X* + 2 C i* I x*|P-l . 
where Bq » 1 / (3 . 


+ X = COS Qt + B 


( 2 . 8 ) 
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2.2.4 CASE IV: Isolation system with asymmetric restoring force 
and subjected to base excitation 

The governing equation of motion can be written as 

M S"+ C 5 + K 2 s^+ 5^= My^ w^cos cat. (2.9) 

Defining the following parameters: 

T = ca^t, n = ca/coQ, A = 

A = S'/ YqCOq, a = S"/y^ ca^ and C = C yg"^ M) , with 

Wq = 

and substituting these in equation (2.9), the following 
nondimensional equation is obtained: 

A + 2 C A |A|P"^ + A + e^A^+ £ 2 ^^"" ' (2.10) 

Where ^ 

By carrying out the transformations along the lines mentioned in 
section 2.2.3, equation (2.10) can be rewritten as 

A* + 2 C A* I A*|^~^ + £2 “ n^cos fix + Bq , (2.11) 

where A* = A + c^/ (3 c^} , B^^ = 1 / (3 and is assumed to 

be satisfy the relation . 
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2,3. Solution Procedure 

Considering the governing equations of motion for the four 
different situations, viz., equations (2.2), (2.4), (2.8) and 

(2.11) , we note that the general form of all these equations is 

Z + 2 C Z i 2 |P"^ + g(Z) = F cos fir + , (2.12) 

where Z, g(Z), F and for the different cases are given in Table 
2.1. In this section, first the procedure for solving equations 

(2.12) in general is outlined. Thereafter, the suitable 
transmissibility expressions for the four different coses (I-IV) 
are obtained. 


TABLE 2.1 


S.No. 

Z 

g(z) 

F 

0 

— 

CASE I 

X 

rH 

1 

X 

X 

\ 

1 

0 

CASE II 

A 


fi^ 

0 

CASE III 

* 

X 

-3 

1 

l/(3v/3r^) 

CASE IV 

A* 

= 2 

fi^ 

l/(3v/3?^) 


The method of harmonic balance 

(2.12). As a first approximation, 
harmonic, and of the form 


[37] is used to solve equation 
the solution is assumed to be 


Z (r) - Zq + 2^ cos (fir + ^) . 


(2.13) 
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Let h(Z,Z) = 2 C Z I Z + g(Z). 


(2.14) 


The periodic function h[Z(T), Z(t)] can be approximated by the 

Fourier expansion up to first order as given below: 

h[Z (r) , Z (r) ] « Rq + cos(nT+0) + R^ sin(Qx+<p) (2.15) 


where the Fourier coefficients R^, R^ and R 2 are given by, 


.Tt 


Rq = 1/(2tt) 


h(il) d-d , 


-TT'' 




R^ = l/( 7T) 


h(il) cos T? di^ , 


(2.16) 


-Tt'’ 




and R 2 = l/( tt) 


-TT' 


h('i>) sin 1 ? dll 


where -d = Qx + <p. 

Substituting equations (2.16) in equations (2.15) and (2.12), 
one obtains 

Z + Rq + R^ cos(nT+(^) + R 2 sin(Qx+(p) = F cos Qr + B^. (2.17) 

Now writing F cos Qr as F [cos{Qx+(f)) cos<p + sin(QT+0) sincj)] 
and equating coefficients of the same harmonics and the constant 
terms from both sides of equations (2.17), one gets 

-Z^Q^ + R^ = F cos Ip 
and R 2 = F sin p, 


(2.18) 
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Where and R„, given by equations (2.16), are functions of 

Zq,Z^, (p and other system parameters. Solving the set of 

simultaneous nonlinear equations (2.18) in Zq, and (p , one 

finally obtains Z(t) from equation (2 . 13 ) . Detailed expressions for 

R , R. and R„ for the four cases (I-IV) are listed in Table 2.2. 
0 1 2 

The symbols Tp and 7^, used in the expressions of R^ and R 2 in 
Table 2.2, are given by 

r = {2/Vn) r([p+2]/2) / r([p+3]/2) (19a) 

hr 

and 

yg = (2/Vn) r([q+2]/2) / r([q+3]/2) (19b) 

Where F is the standard gamma function. 

2.3.1 Transmissibility expressions 

A minor modification is made in the standard definition of 
the transmissibility to account for the constant term present in 
the response of the system with asymmetric restoring force. Two 
indices, namely the force transmissibity (T^) in the case of force 
excitation and the displacement transmissibility (T^) in the case 
of base excitation, are defined as given below. 

The force transmissibility is defined as the ratio of the 
maximum magnitude of the force transmitted to the foundation to 
that of the exciting force (F) . From equation (2.17) it can be 
seen that the force transmitted to the foundation is given 
by, 


Ft = Rq + cos(Qx+<P) + sin(Qz+<P) - b^. (2.20) 


Thus, we get 
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TABLE 2.2 
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Th6 displacGinent transmissibility dofiticd as the 

ratio' of ths inaxiHiuin displacsinont of M to that of the bnbo iinotion. 
Aftor al^sbraic manipulation/ it can bo shown that tho ^onBiralized 
expression for is given by 

= ![ Zq - 3Bq + (zJ + 1 + 2 Z^ cos(p ]|. (2.22) 

Detailed expressions of and for cases I to IV are given 
below. These quantities can be evaluated by referring to Table 

2 . 2 . 

Case I: = [K^ + 

[(K^ - + cl ^2.23) 

where = 2 C 12^“^ y 

e ^ 1 "p 

and . xf"- . , 2 . 24 , 

It may be noted that the method of equivalent linearisation 
[22,110] gives the same results as obtained above, where, C and 
are referred to as the equivalent viscous damping and the 
equivalent stiffness, respectively. 

Case II; is still given by equation (2.23), but in equation 

(2*24) should be replaced by 


case III: = ICR^ - 1/(3 + (r2 + r2,1/2 

^ 1 2 


(2.25) 
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Case IV; 


[ - 




(A^ + 1 


+ 2 


A* COS0 ) 


]l 

(2 .26) 


2.4 Linear Stability Analysis 

In equation (2.12), both restoring and damping forces do not 
have continuous derivatives with respect to Z and Z, respectively, 
for p < 1 or/and q < 1. Hence, the stability analysis by the usual 
variational technique is not amenable for analytical treatment. A 
modified procedure described in reference [111] is, therefore, 
adopted. Considering that Z^ and Z^ in equation (2.13) are slowly 
varying with time and substituting it in equation (2.17), the 
assosiated autonomous equations of the system are obtained as 
given below: 

i - Fcos 0 ) - Z^ 

f 

2 z^ n 

Z^ = -(Fsin <p - )/{ 2 Q ) , (2.27) 

and = -Rq “Bq 

While deriving equations (2.27), Z^ but not Z^ has been neglected 
[42]. The relevant expressions for Rg , and R 2 can be 

substituted from Table 2.2 for different cases I-IV. The stability 
of the solution, 0=:Z^ = Zg=^ = O can be studied by considering 
the linearised equations 


7} = A T) 


(2.28) 
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Where t] = { <t>, Z^, Z^, ^ and A is the Jacobian matrix of the 
system under consideration. The eigenvalues of A provide the 
requisite information on the stability of the solutions obtained 
in section 2.3. When both p and q are odd integers, the present 
method and the usual variational technique yield the same results. 
An alternative approach to carry out the stability analysis using 
the notion of "equivalent viscous damping" is proposed in the 
appendix A for the case of Buffing's oscillator with pth-power 
damping. 

2.5, Results and Discussions 

In this section a parametric study depicting the effects of 
system parameters on the performance indices of various types of 
isolators is presented. Four typical values of p = 1, 1.5, 2, 3 
signifying, respectively, viscous, orifice, quadratic and cubic 
damping are considered. Results for a symmetric restoring force 
are reported in section 2.5.1 and those for an asymmetric 
restoring force are presented in section 2.5.2. A representative 
curve indicating transmissibility ( in dB ) vs frequency (Q) is 
shown in Figure 2.3. 

With refejfrence to this figure, the following indices are 
defined : (i) the maximum value of transmissibility (T^) , (ii) the 
frequency at which the maximum transmissibility occurs (Q ), (iii) 
the resonance critical frequency (Q^) , (iv) the non-resonance 
critical frequency (n^) , (v) the jump width (W) and (vi) the high 
frequency attenuation rate (HAR) of transmissibility. The 
existence of three possible values of the transmissibility (Figure 
2.3) at a given frequency gives rise to the so-called jump 
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phenomenon. The upper branch of the jump is known as the resonance 
the lower one as the non-resonance branch and the middle 
one as the unstable branch. Two points of vertical tangoncies on 
the transmissibility curve are referred to as the critical 
freguencies. One on the resonance branch is identified as 0^, 
while that on the non-resonance branch is indicated as Q^. The 
jump width (W) is given by W = . An unstable solution of 

the governing equation exists in the frequency range given by the 
jump width. The high frequency attenuation rate of 

transmissibility is given by the slope of the transmissibility (in 
dB ) vs curve (as Q tends to oo) . 

It can be seen from Figure 2.3 that the difference between 

and is negligible, and one can consider and the 

associated transmissibility (T^) on the resonance branch as n and 

r m 

T^, respectively. Obviously a good isolator should have low values 
of and T^ and a high value of HAR. 

2.5.1 Isolation system with symmetric restoring force 

The transmissibility expression for a force excitation is 
given by equation (2.22), and that for a base excitation by 
equation (2.24). The transmissibility values of a soft system 
(q<l) for these two excitations are shown, respectively, in 
Figures 2.4 and 2.5. It can be seen from these figures, that a 
jump similar to the one observed in the response of a soft 
Duffing's oscillator, occurs. However, a soft Buffing's 
oscillator, in the absence of excitation, possesses three 
equilibrium points [36] whereas the present system in the unforced 
has only one equilibrium point. Therefore, an anomalous jump 
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observed in the response of a soft Duffing' s oscillator with high 
excitation values [74] is absent in the present system. For a hard 
system (i.e., q>l) , a jump similar to the one observed in the 
response of a hard Duffing' s oscillator occurs as shown in Figures 
2.6 and 2.7. In general, an isolator is effective if the resonance 
frequency ( £1^ in this case) is far below the operating range of 
interest. One may observe (see Figures 2. 6-2. 7) that the jump 

width for a hard system covers the high-frequency range, and 

consequently such an isolator is not as effective as an isolator 
with a soft characteristic. 

The effect of increasing ^ can be summarised (See Figures 

2. 4-2. 7) as follows: With increasing C/' (i) the resonance 

transmissibility decreases, (ii) the high frequency attenuation 

rate (HAR) decreases, (iii) the jump width (W) reduces and may 

eventually be entirely eliminated, (iv) the resonance frequency 

(n^) decreases for q > 1 and increases for q < 1. Thus, the 

effects of C on T & HAR are similar to those for a linear 

r* 

isolator. 

The jump width (W) reduces as the damping exponent p 
increases as can be seen from Figures 2.8 and 2.9. It is 

interesting to note that is insensitive to variations in p and 
Thus, in many cases one can arrive at an analytical estimate of 
by approximating C = 0* As expected, the resonance critical 
frequency is very much sensitive to variations in the damping 
parameters p and <. Furthermore, the jump width in the case of a 
base excitation is more than that for a force excitation when q > 
1 and vice versa when q < 1. 

Referring back to Figures 2.5 and 2.7, one sees that HAR 
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decreases with an increase in p for the case of a base excitation. 
It becomes zero for p=2 and even negative for higher values of p. 
Comparing Figures 2.4 and 2.6 with Figures 2.5 and 2.7, 
respectively, it is seen that the effects of p on HAR for a force 
excitation is just the opposite of that for a base excitation. 
With a force excitation, HAR increases with increasing p. Thus, 
the trends of HAR are similar to those reported by Ruzicka [22] 
for the case of an isolator with linear stiffness. 

The stability analysis is carried out as explained in section 
2.4. The unstable branch in the transmissibility plot is shown by 
dotted lines in Figures 2. 4-2. 7. The eigenvalues of the Jacobian 
matrix A given by equation (2.17) show that this unstable solution 
corresponds to an unstable saddle point of the associated 
autonomous equation given by equations (2.16). 


2.5.2 Isolation system with asymmetric restoring force 

The transmissibilities given by equations (2.25) and (2.26) 
are plotted in Figures 2.10-2.11. Numerical computations are 
carried out with = l.o, C, = o.Ol and 0.5,. A jump similar to 
that in a 'hard' Duffing's case is observed in all these figures. 
The unstable branches are marked with dotted lines. The effects of 
^ and p on all the quantities associated with the jump are similar 
to those observed with a symmetric restoring force. But the 
effects of C and p on HAR with a base excitation are at variance 
from that of the previous section. In the present case, HAR always 
tends to zero. For the parameter values considered here, one can 
conclude that an isolation system with asymmetric restoring force 
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may not perforin satisfactorily with a base excitation. 

2.6. Conclusions 

An isolator with a soft characteristic (q<1.0) is found to 
be superior to that with a hard characteristic (q>1.0), since in 
the latter case the jump in the transmissibility curve extends 
towards high frequencies. The effects of damping on the resonance 
transmissibility and the high frequency attenuation rate are 
similar to those reported by Ruzicka for systems with a linear 
restoring force. The jump width decreases with increasing damping 
coefficient and may even entirely be eliminated. With increasing 
damping index p, the jump width reduces making the unstable zone 
narrower. The resonance critical frequency is very much sensitive 
to variations in the damping parameters whereas a good 
approximation for the non-resonance critical frequency can be 
obtained by assuming the damping to be absent. It is seen that for 
the parameter values considered in the present work, an isolator 
with an asymmetric restoring force may not perform satisfactorily 
for a base excitation. 



CHAPTER 3 


CHAOTIC RESPONSE OF NONLINEAR ISOLATORS WITH pth-POWER DAMPING 

3.1 Introduction 

The steady-state, harmonic response and the transmissibility 
indices of nonlinear isolators have been discussed in chapter 2. 
In the present chapter, chaotic oscillations of a harmonically 
excited mass on a nonlinear isolator are investigated. The 
magnitude of the strictly dissipative damping force on the mass, 
provided by the isolator, is assumed to be proportional to 
pth-power of velocity. Both symmetric and asymmetric nonlinear 
restoring forces are considered. Numerical simulation is carried 
out to analyse the subharmonic and chaotic motions in the 
background of the results obtained in the case of linear damping 
(i.e., p=l) . Two typical routes to chaos, namely through the 
period-doubling and intermittency , reported already in the 
literature for linear viscous damping (i.e. p=l) , are seen to be 
present with the damping exponent p=2 and p=!3. Thus, the 
bifurcation structure seems to be unaffected by the damping 
exponent p. Of course, the values of the damping coefficient 
(constant of proportionality) needed for complete elimination of 
the subharmonic and chaotic responses depend on the value of p. A 
parametric study is presented to indicate the role of the damping 
exponent, damping coefficient and asymmetry on the onset of chaos. 
It is revealed that suitable choices of the damping coefficient 
and damping exponent can entirely eliminate subharmonics and 
chaotic motions and thus provide a passive rather than an active 


control scheme for quenching chaos. 
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3.2. Isolator Model 

The governing equation of motion, for a harmonically excited 
mass on an isolator with pth-power damping and cubic stiffness 
nonlinearity, can be written as 

X" + 2Cx'lx'|^“^ + x^ = f^ + f cos wt; C>0/ P>0 (3.1) 

where x is the displacement of the mass and o) is the excitation 
frequency. It may be noted that equation (3.1) has been obtained 
by following a non-dimensional procedure different from that used 
in chapter 2. In this chapter, results are obtained by varying 
excitation with cj=l.O. The values of the damping exponent p = 1, 2 
and 3 correspond, respectively, to linear, quadratic and cubic 
damping models. The case p=0 representing Coulomb damping with a 
discontinuous damping force has many special features and will be 
presented in a later chapter. The parameters ^ and f in equation 
(3.1) refer to the damping coefficient per unit mass, and the 
amplitude of the exciting force per unit mass, respectively. The 
constant term f^ represents the asymmetry of either the excitation 
or the restoring force. In a real-life vibration isolation system, 
this asymmetry can arise due to various factors such as when (i) 
the gravity effects are included in a vertical system (ii) the 
restoring force has a quadratic term as in a constant area 
diaphragm pneumatic spring. In the latter, a suitable scaling [37] 
gives rise to the constant term f^ in the excitation- Overall, the 
degree of asymmetry in the system is represented by which plays 

central library 

> »• T,, KANP UH 

A. J 
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a major role in the analysis of chaotic motion as will be pointed 
out later. In this chapter, the symmetric (1^=0) and asymmetric 
cases are discussed separately. It should be noted that a 
linear term which may be present in the restoring force of a 
typical isolator has been omitted in equation (3.1) since it is 
known that the linear term does not qualitatively affect the 
behaviour of the system in the chaotic regime [53-55]. 

As indicated by the enormous complexity of the results 
obtained by numerical simulation over the last ten years [43-71], 
it is a hopeless task to compute the transmissibility indices for 
all the parameter zones of interest. Consequently, in this chapter 
the response of the mass rather than any isolator index has been 
computed to obtain the qualitative behaviour of the system, 

3.3 Numerical Results and Discussions 

In this section, detailed numerical simulation of the system 
governed by equation (3.1) is presented. The effect of nonlinear 
damping on the bifurcation structure and routes to chaos are 
discussed. Towards this end, the usual phase plots, Poincare maps 
and response plots obtained from numerical simulation are used. 
The Runge-Kutta-Merson method has been employed to carry out the 
numerical integration of equation (3.1). 

t 

3.3.1 Symmetry-breaking and period-doubling route to chaos with 
symmetric nonlinearity 

It should be noted that equation (3.1) has the symmetry such 
that X > -X as t — > t+n/o). In other words, if x(t) is a 
solution of equation (3.1), then so is -x(t+TT/u)) . These solutions 
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may or may not be distinct. The solution is called symmetric if 

X(t) = -X (t+7T/0j) . (3.2) 

It is obvious that any solution consisting of only odd 
(super) harmonics is symmetric. It can be shown [56, 59] through 
the stability analysis (using Mathieu's equation) that at certain 
critical parameter values, the symmetric solution becomes unstable 
giving rise to even order superharmonics. It may be noted that as 
soon as even superharmonics appear, the response loses symmetry, 
i.e., x(t) and -x(t+Tr) solutions become distinct. This loss of 
symmetry manifests itself in the appearance of dual solutions [62, 
67] as confirmed by numerical integration with different initial 
conditions. With further increase in f, this period-doubling route 
(i.e., through the introduction of cj/2 , w/4, ... in the response) 
eventually leads to chaotic response. It has been shown in the 
literature that the symmetry-breaking precursor is necessary for 
the period-doubling route to chaos in symmetric systems [61]. 

Numerical simulation of equation (3.1) is carried out with 
the following set of parameter values: u- 1.0, 0.025, p=2 . 0 and 
fQ=0. Bifurcation sequences are observed by varying the amplitude 
of excitation, f. 

The phase plots, obtained after the transients have died 
down, are shown in Figure 3.1. In Figure 3.1a, the time period of 
the periodic solution is found to be the same as that of the 
excitation. The symmetry of the phase trajectories suggest that 
only odd superharmonics appear in the solution. In Figure 3.1b, 
the symmetry of the phase trajectory is lost and the dual solution 
(for the same values of all the parameters) obtained with a 
different set of initial conditions is shown in Figure 3.1c. Both 
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these solutions have time period equal to that of the excitation. 
Figure 3 . Id shows the phase trajectory when 1/2 subharmonic 
appears in the solution and consequently the period of the 
response is twice that of the excitation. Figures 3.2a and 3.2b, 
respectively, show the Poincare map and the response plot of the 
chaotic attractor obtained through this period-doubling route. 

3. 3. 1.1 Effect of damping exponent p and damping coefficient C 

A parametric study is carried out to determine the 
effects of p and ^ on the onset of chaos via the period-doubling 
route. Towards this end, the critical forcing amplitudes are 
obtained for the occurence of symmetry-breaking and the first two 
period-doubling bifurcations. Figures 3.3a-3.3c refer to the 
results obtained with p=l, 2 and 3, respectively. It can be seen 
from Figure 3.3a that a minimum value of C is necessary for 
breaking the symmetry. Further, for a given value of p, the 
forcing amplitude required to break the symmetry is relatively 
independent of the value of C* However, as C, increases there is a 
marked increase in the forcing required for the successive 
period-doubling bifurcations. It may be noted that with increasing 
C, the cascade of period-doubling bifurcations reduces and chaos 
is entirely suppressed though the initial period-doubling may 
still continue. Hence, strictly dissipative nonlinear damping, 
just like linear damping, can also be used as a passive mechanism 
to suppress chaos. 

Comparing Figures 3.3a-3.3c one can conclude that the basic 
bifurcation structure is totally independent of the value of p. 
Only the minimum value of needed for the symmetry-breaking, or 
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Figure 3,3 Effect of p and C on the critical forcing amplitude 
required for symmetry-breaking and period-doubling (a) p=1.0; (b) 
p-2,0; (c) p=3,0; Symmetry-breaking Period 2 ; Period 4 






58 


^Qj- ths suppression of period-doubling cascddes (and eventual 
route to chaos) decreases with increasing p. The designer of 
isolation systems can suitably choose p and C to eliminate chaos. 

3.3.2 Intermittency route to chaos with symmetric nonlinearity 

In the present section, another route to chaos without 
period-doubling known as intermittency will be illustrated for 
p=2, through the numerical results obtained by integrating 
equation (3.1). This chaotic response is seen to be associated 
with a transition from a period 3 orbit to a period 1 orbit. 
Detailed mathematical investigations of chaos arising out of 
period 3 solution has been presented by Li and Yorke for 
one-dimensional maps [117]. The intermittency transition has been 
introduced by Pomeau et. al based on the numerical simulation of 
Lorenz system [118]. Since then, numerous papers have appeared 
introducing various types of intermittency transitions to chaos 
[43-46, 119]. 

The transition from a period 3 orbit to a period 1 orbit via 

a chaotic zone has been observed in a Puffing's oscillator with 

linear damping [53, 70] and has been called the sharp transition 

to chaos in reference [70] , where a detailed model is given based 

on the 1/3 and 7/3 subharmonics. Both Type I and Type III 
• • ^ 

intermittency routes have been observed in Duffing-type 
oscillators [68]. Detailed literature on the renormalisation 
procedures in intermittency is available and incidentally this 
happens to be the only route where exact renormalisation is 
possible [46]. Even in the logistic map, beyond the critical value 
r^ (of the Feigenbaum period-doubling route) , period three windows 
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occur which either go to chaos through period-doubling or through 
Type I intermittency [45] . It is interesting to point out that 
such a period three orbit and subsequent period-doubling route 
have also been observed in Duffing's equation [62, 67]. The 
numerical results of driven dissipative oscillators thus resemble 
the logistic map more readily than the circle map [64], 

Using the same values of the parameter given in section 3.3.1 
(i.e., <j)= 1.0, 0.025, p=2.0) and increasing the value of f, a 
period 3 orbit is revealed whose Poincare map is shown in Figure 
3.4a. With further increase in f, a prechaotic zone is obtained. 
The Poincare map of this zone consisting of three small line 
segments centered around the period 3 orbit is shown in Figure 
3.4b. The Poincare map of the chaotic solution obtained with 
further increase in f is shown in Figure 3.4c. It can be seen from 
Figure 3.4d that further increase in the forcing results in a 
period 1 orbit. The chaotic response is essentially a transition 
from a period 3 orbit to a period 1 orbit. Thus, even the 
intermittency route to chaos also appears to be independent of the 
value of p. 

3.3.3 Asymmetric isolation system 

As mentioned earlier, the asymmetry of the isolation system 
reflected by the constant forcing term fQ in equation (3.1) plays 
a major role on the onset of chaos. Results of section 3.3.1 
indicate that for a symmetric system, the breaking of symmetry is 
essential for period-doubling. Further, the symmetry-breaking 
yields a dual solution and both the solutions go to chaos, through 
period-doubling, resulting in dual chaotic attractors [67]. When 
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asymmetry is inherent in the system, the symmetry is trivially 
broken or rather the symmetry-breaking bifurcation is not 
necessary for period-doubling. Dual solutions appear right from 
the beginning. Both the solutions undergo period-doubling 
bifurcations. But one of these routes gets excited for a very 
small value of the forcing amplitude f even for a slight degree of 
asymmetry (fQ) • 

The above scenario is seen to occur for p=l [54]. In the 
present section it is shown that similar results are obtained even 
for p=2. Numerical simulation of equation (3.1) is carried out 
with the following parameter values: fQ= 0.03, cj= 1.0, 0.025, 

p=2.0. The Poincare map of the resulting period 2 orbit is shown 
in Figure 3.5. With increasing values of f, a chaotic attractor is 
obtained (via the period-doubling route) whose Poincare map and 
time response are shown, respectively, in Figures 3.6a and 3.6b. 
Comparision of the values of f required for the onset of the 
period 2 orbit in Figures 3 . Id and 3.5 confirms that a little 
asymmetry drastically reduces the necessary forcing amplitude for 
beginning of one of the period-doubling routes to chaos. 
Similarly, from Figures 3.2b and 3.6b, it may be noted that the 
amplitude of the resulting chaotic attractor along this route is 
small. However, in this asymmetric case, there exists another 
period-doubling route to chaos for a different set of initial 
conditions. This route, shown in Figures 3.7a-3.7c, is relatively 
unaffected by the degree of asymmetry in the system. Along this 
route, the value of f required to cause the period-doubling and 
the amplitude of the chaotic attractor are more or less of the 
same magnitudes as those of the symmetric case. The intermittency 





80*0 = 













65 


route to chaos, explained in section 3.3.2 for the symmetric case, 
also occurs in the asymmetric case and is shown in Figures 
3 .83-3 . 8d. 


3.3.4 Discussion 

All the above scenarios presented for quadratic damping are 
found to exist irrespective of the value of the damping exponent 
p. Even combinations of various types of damping viz., quadratic 
and linear or linear and cubic etc., yield similar trends. It is 
seen that the bifurcation structure associated with the 
resonances of the system [62-69] is not changed by such nonlinear 
damping models so long as they are strictly dissipative. However, 
as reflected by the parametric study, the critical values of 
forcing at which various bifurcations occur depend on the values 
of p and C* Thus, a proper selection of the damping mechanism with 
suitable parameters provide the designer with wider choices for 
passive control of chaos. The fact that the bifurcation structure 
is not altered by the value of p also points to the possibility of 
extending the notion of equivalent linear damping to obtain the 
approximate criteria of chaos based on the stability of harmonic 
solution. 

The results of only the primary resonance region are 
presented in this work. It is known that the secondary as well as 
subharmonic resonances exist in these systems and the 
period-doubling and intermittency routes recursively repeat 
themselves [62-70]. The results presented in this chapter are for 
a fixed frequency (w=1.0) with varying values of f. However, as 
found earlier with linear damping [70], similar trends are 
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observed for a fixed value of f with varying frequency even with 
nonlinear damping models. 

3. 4 Conclusions 

Vibration isolators with nonlinearity in stiffness and 
damping terms are analysed under harmonic excitations. It is shown 
that both the period-doubling and intermittency routes to chaos 
occur for all values of p, thus preserving the bifurcation 
structure of a linearly damped, buffing's oscillator. A parametric 
study is reported to indicate the effect of damping on the 
symmetry-breaking and period-doubling. It is pointed out that 
asymmetry inherent or arising out of gravity effects in an 
isolation system plays a major role on the onset of chaos. It is 
shown that nonlinear damping, just like linear damping, can also 
be effectively used for controlling chaos. 



CHAPTER 4 


HARMONIC RESPONSE OF HARD DUFFING-TYPE VIBRATION ISOLATOR WITH 
COMBINED COULOMB AND VISCOUS DAMPING 

4.1 Introduction 

In the present chapter, the steady-state, harmonic response 
of a vibration isolation system with a cubic, hard nonlinear 
restoring force and combined Coulomb and viscous damping is 
presented. The results have been obtained by using the method of 
harmonic balance. It has been assumed that the motion is 
continuous without any stop. Both force and base excitation models 
are presented. 

While analysing the steady-state response of a viscously 
damped, soft Duffing-type system by using the method of harmonic 
balance, an anomalous jump in addition to the standard one has 
been observed [72-80] . The importance of this anomalous jump in 
understanding the so-called symmetry-breaking phenomenon in soft 
systems has been discussed in references [75-77]. It is 
interesting to note that a similar anomalous jump has been 
observed in the present work while analysing the steady-state 
response of a base excited hard, Duffing-type system with Coulomb 
damping. A linear stability analysis following the method 
suggested in reference [111] is carried out to ascertain the 
status of this additional jump response. Also some analytical 
criteria for merger of the normal jump with the anomalous one are 
derived. The effect of dissipation on this anomalous jump has been 
studied in detail. The existence of an anomalous jump in the 
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response obtained by the method of harmonic balance, is confirmed 
by direct numerical integration. However, unlike in the case of a 
soft system, this anomalous jump seems to have no bearing on the 
symmetry-breaking and subsequent period-doubling route to -chaos. 

The transmissibility indices have been plotted for relevant 
parameter values to study the performance characteristics. 
However, additional features like the effects of negative damping 
and response with stops have not been addressed to in this thesis. 
The results of this work extend those obtained previously by Den 
Hartog and Ruzicka for a linear spring and may thus enhance the 
scope of the design of vibration isolation systems. 

4.2 Theoretical Analysis 

4.2.1 Isolation system with base excitation 

Referring to Figure 4.1, let a rigid mass M be separated from 
a harmonically moving base through an isolator consisting of a 
symmetric, nonlinear spring and combined Coulomb and viscous 
damper. The equation of motion for the mass can be written as 

M 5" + C 5' + sgn(6' ) + K 5^ = -M y" = My^ cos wt 

( 4 . 1 ) 

Where C and are the viscous and Coulomb damping coefficients, 
respectively; 5 = (x-y) is the relative displacement of the body 
to be isolated, and the prime denotes differentiation with respect 


to time t. 
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Figure 4.1 Vibration isolation system with combined Coulomb and 
viscous damper and a cubic nonlinear spring subjected to base 
excitation 



viKous'^dLoIr^at^^^" isolation system with combined Coulomb and 
IxKtation ® nonlinear spring subjected to force 





71 


Defining the following non-dimensional parameters: 

T = Wq t, n = W /Wq, 

A = S/Yq , A = S'/CYo ^ = S'VCYo )/ C = C/(2MWq), 

Cf = C^/(2MyQWQ2) with = [ K and 

substituting them in equation (4.1), one obtains, 

A + 2CA+2Cf sgn(A) + A^ = cos Qz (4.2) 

where the dot denotes differentiation with respect to 
non-dimensional time x. 

Assuming the steady state solution of equation (4.2) to be given 
by A(x) = A^ cos (Qz+(p) and using the method of harmonic balance, 
the following equations are obtained: 



-A^Q^ + 

2 

= Q cos 

0 

(4.3a) 

and 


sin 0 


(4.3b) 

where 


Of 


(4.4a) 


^1= 

-(2CA^n + 2 

Cf 13 ) 

(4.4b) 


with a- 3/4 and 4/ti 


Eliminating <p from equations (4.3a) and (4.3b), one can write the 
amplitude equation as 


A^^ -2 a Q^A^'^+ (Q^ + 4 + 80C<^nA^ 

+ (4/3^ Cf^- = 0- 


(4.5) 
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In the absence of viscous damping (i.e.,with C-0) / equation (4.5) 
reduces to 

-2 a A^^^ + (4/3^ (:/ - Q^) = 0. (4.6) 

Some simple analytical criteria are derived in the Appendix C for 
the existence of repeated roots of equation (4.6). The 
significance of the repeated roots of equation (4.6) will be clear 
later on. 

One can observe that if the constant term in equation (4.5) 
(same as that in equation (4.6)) is equated to zero, then one root 
of A^ becomes zero implying no relative movement. Thus, the 
well-known break-loose frequency (Qj^) can be obtained by equating 
this constant term to zero as: 

4 13^ c/ - = 0 (4.7) 

i.e, = (2 0 . (4.8) 

When the excitation frequency is below the break-loose frquency, 
the damper gets locked and the relative motion does not take place 
(i.e., 5 H 0) . Hence, the present approximate analysis based on 
continuous relative motion is not valid below the break— loose 
frequency. 

Finally, in the absence of Coulomb damping (i.e., with Cf “0)/ 
equation (4.5) reduces to 

- 2 a A^^ + (fi^ + 4 0^) A^^ - = 


0 . 


(4.9) 
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. . . . . 2 
Equation (4.9) is a cubic equation in which is nothing but the 

well-known solution of a viscously damped buffing's oscillator.lt 

may be noted, that the method of equivalent linearization [110] 

gives the same results as obtained above with the equivalent 

stiffness and viscous damping parameters defined by 

K = a A-^ 
e 1 

and = 2C + 2Cf^/(A^fi). (4.10) 

The absolute displacement transmissibility (T^) is defined as 
the ratio of the maximum displacement of the mass to that of the 
base motion. It can be shown that 

= [1 + A^^ + 2 A^ cos (4.11) 

where A^ is given by equation (4.5). Substituting for cos <p from 
equation (4.3a) in equation (4.11) one can obtain 

= [ (n^ + 2aA^'^ - A^^ )/(f2^) (4.12) 

Alternatively, following the linear vibration theory, one can 

t 

write [3] 

T, = [K ^ + C J / [(K -Q^)^ + C ^ (4.13) 

Substitution for K and C from equations (4.10) in the above 
equation results in 
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Td = [t^ / t2 (4*14) 

where = oi} + 4 

and t^ = oc^ -2 a n^A^'^+ (Q^ + 4 + B^CCf^A^ 

2 2 

+ Cf - 

with the help of equation (4.5), it is easy to verify that 
equations (4.14) and (4.12) are identical. 


4. 2. 1.1 Stability analysis 

It can be seen from equation (4.1) that the derivative of the 
damping force with respect to S' does not exist at 6' = 0. Hence, 
in place of the usual variational technique, a modified procedure 
for the stability analysis as suggested in reference [111] is 
followed. Using this approach, the stability of the steady state 
response is predicted by the stability of the equilibrium points 
of the following set of associated autonomous equations: 

(p = (aA^^ - cos < 1 ^ “ A^ n^)/(2A^n) 

(4.15) 

* 2 

A^ = -(Q sin (f) + 2CA^Q + 2C^^)/(2n). 


The stability of the equilibrium points of equations (4.15) 
can be studied by considering the equations 


7} = J 7) 

where 7}={«;J, A^ }'^ 


(4.16) 

and J is the Jacobian matrix whose eigen values. 
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when evaluated at the equilibrium points of equations (4.15), 
provide the requisite information on the stability of the 
solutions obtained [42]. It is easy to see that the matrix J is 
given by 


J = 




(4.17) 


where = n sin (/>/(2A^), 

a A^/n + n cos 0/(2 A^) , 

cos 0 /2 (4.18) 

and *^4 • 


4.2.2 Isolation system with force excitation 

Referring to Figure 4.2, when the rigid mass M, acted upon by 
a harmonic force, is separated by the isolator from an immovable 
foundation, its equation of motion can be written as 

M x' ' + C X' + Sgn(x') + Kx^ = f cos wt. (4.19) 

Defining the following non-dimensional parameters: 

T = Wq t, n = W/Wq, 

X = X/Xq , X = xV(XqC0q ), X = x'VCXq ^ C/(2MWq), 

Cf = C^/(2MxQWg^) with Xq = (f/K)^/^, Wq = [K Xq^ ]V2y j^l/2 

and using them in equation (4.19), one obtains, 

X + 2CX+2Cf sgn(X) + X^ = cos fJx . (4.20) 

Assuming the steady state harmonic response to be given by X(t) = 
X^ cos (nT+0) , the response equations for X can be obtained as 
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(4.21a) 

and 

= sin (p 

(4.21b) 

where 

R2= d 

( 4 . 22a) 


S^= "{2CX^n + 2 Cf /3) 

(4.22b) 

with 

a= 3/4 and (3= 4/tt . 



% 


Now, using equations (4.21) and (4.22), the amplitude equation is 
finally obtained as 

oi^ -2 a (n"^ + 4 + S/lCCffiX^ 

+ ( 4^2 = 0 ^ ^^ 23 ) 

In the absence of viscous damping (i.e.,with C=0) , equation (4.23) 
reduces to 

Xi® -2 a X^2 + (4^2 Cf^ - 1) = 0 . (4.24) 

Denoting the critical value of at which the constant term in 

equations (4.23) and (4.24) goes to zero as (Cj)^, one can see 
that 


*^f^o ■ i -0 (4.25) 
or, (Cf)^= 1/(2^) = r,/8 . 

Kith <=(Cj)^ one root of the amplitude equation goes to zero. As 
Shown )5Y Den Hartog [91] for the isolator with a linear spring and 
a coulomb damper, equations (4.23) and (4.24) are not valid if C, 
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> • Thus, the approximate solution obtained by the method of 

harmonic balance is valid for Cf < tx/S. The force transmissibility 
(T^) is defined as the ratio of the amplitude of the force 
transmitted to the foundation to that of the exciting force and is 
given by 

Tf = /f = [t^/t^]^/^ (4.27) 

Where t^ = + 4 + 8/3CCffJX^ + 4/3^ 

and X^® -2 a + 4 + 8/3CCfS^X^ 

+ 4/3 Cf • 

Comparing equations (4.14) and (4.27) one can conclude that 
equation (4.27) is nothing but equation (4.14), when X^ is 
substituted for A^. The stability analysis of the steady state 
response is carried out in the same manner as followed for the 
base excitation model. 

4.3 Results and Discussion 

Performance characteristics of the isolation systems 
described in section 4.2 are presented through the 
transmissibility in dB (i.e., 20 log^^ or 20 log^Q T^ as the 
case may be ) versus frequency (fi) plots for various values of ^ 
and Cf 

4.3.1 Base excitation 

Case(i) In the absence of viscous damping (i.e. (^=0) : 

The absolute displacement transmissibility (T^) is given by 
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equation (4.14) where is evaluated from equation (4,6). a 

typical transmissibil ity versus frequency (T^ vs Q) plot is shown 
in Figure 4.3. One can define all the parameters relevant to the 
present discussion with respect to this figure. As mentioned 
earlier, when the excitation frequency is below the break-loose 
frequency (ilj^) , the Coulomb damper gets locked and the 
transmissibility is maintained at unity. This portion is denoted 
as OA in Figure 4.3, with point A corresponding to the break-loose 
frequency The resonance branch of the usual Duffing jump is 

labelled DBF, the unstable branch as IG and the non-resonance 
branch as IH. The frequency corresponding to the point I is 
denoted as non-resonance critical frequency (Q^) . It can be seen 
from Figure 4 . 3 that the stable resonance branch DBF extends 
towards infinity resembling the response of an undamped 
oscillator. 

In the case of a viscously damped, soft Duffing oscillator 

with a high value of excitation, an anomalous response apart from 

the usual jump phenomenon has been reported [72-80]. It may be 
noted from Figure 4.3 that a similar peculiar behaviour is 
observed for the system under consideration when Q > Just 

beyond the break-loose frequency, a small portion with two 

additional branches (labelled BC and AC in Figure 4.3) appear 
along with the resonance branch DBF. This portion with two 

branches BC and AC is denoted as the 'anomalous jump' . Let the 
frequency at which these two branches meet (at C in Figure 4.3) be 
denoted by The range of frequency over which this anomalous 

jump exists depends on the value of 

At a given value of Cf/ equation (4.6) has repeated roots 
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When n = n and n = n . It is of interest to examine the effect of 

C on the values of Q and Q^, A procedure for obtaining the 

condition for repeated roots of equation (4.6) is explained in the 

Appendix C. Solving equation (C.6) one can obtain and for 

various values of Cf These results along with the variation of 

with respect to Cf (refer to equation (4.8)) are plotted in Figure 

% 

4.4. Upon increasing C,^, the range of frequency (Q^-Qj^) over which 
the anomalous jump exists also increases and decreases. At 

a critical damping value Cf=(Cf)jn ( indicated by the point M in 
Figure 4.4) equals and the anomalous jump coincides with the 
main jump. At this critical value (Cf)j^ equation (C.6) has 
repeated roots . The expression for (Cf)j^ is also derived in the 
Appendix C. 

The nature of variation of the transmissibility with 
frequency as Cj increases is shown in Figures 4 . 5 (a) -4 . 5 (d) . It 
may be noted from Figures 4.5(a) and 4.5(b) that with low values 
of the anomalous portion is very small. This is not surprising 
since the difference Qj^) is very small with low values of Cf 

as can be verified from Figure 4.4. For Cf > (Cf)jn' 
branch of the usual jump opens up and merges with the anomalous 
jump as clearly seen in Figure 4.5(d). 

Now referring again to Figure 4,3, with a low value of Cf < 
as the frequency is increased beyond the response takes 
the lower branch AC of the anomalous portion and at Q=n (i.e., at 
the point C) it jumps to the point E on the usual resonance branch 
DEF • On the otherhand, if the isolator is operating on the 
non-resonance branch IH, then with decreasing frequency the jump 
takes place from the non^-resonance branch to the resonance branch 



Figure 4.4 Variations of Q^, and with in the absence of 

viscous damping (C==0) . 0 ; Q q . 
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at i.e., from the point I to F, and the response follows 

the curve FED. The response again jumps back to the locked region 
at the point D. 

It is known that in the case of a linear spring and a 
Coulomb damper [3], with a suitable choice of the break-loose 
frequency, (i.e., with a friction coefficient > n/4) , one can 

maintain the transmissibility less than or equal to unity. Also 
the unbounded transmissibility at resonance can be avoided if the 
break-loose frequency is made atleast equal to 1 ( i.e, with a 

friction coefficient > tt/S ). In the present case, by choosing Cf 
> shown in Figure 4.5d, one can maintain the 

transmissibility to be equal to or less than unity over all 
frequencies as no jump occurs towards the resonant branch. It can 
be seen from the Appendix C ( refer equations C.7 and C.8) that 
7T/4 and the corresponding break-loose frequency • 

The stability analysis of the anomalous jump portion (see 
Figure 4.3) is carried out in a manner detailed in section 4.2.2. 
It is found out that near the break-loose frequency with Qj^< Q < 
the. upper and lower branches i.e. AC and DE are either stable 

a 

nodes or stable spirals and the middle branch (i.e.,BC) 
corresponds to unstable saddle points. Due to the transition from 
a stable node to a stable spiral near Q=Qj^, some interesting 
transient behaviours are expected. 

Case(ii) Effect of viscous damping: 

The effect of on the anomalous jump in the presence of 
viscous damping can be seen in Figures 4. 6-4. 8. The break-loose 
frequency is of course independent of ^ (see equation (4.8)). 
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choices of C Cf made for the desired performance of 

the isolation system according to the application on hand. Results 
obtained, not presented here, with an asymmetric restoring force 
and Coulomb damping also show the anomalous jump. 

4. 3. 1.1 Numerical integration results 

In the present section, equation (4.2) is integrated 
numerically (using Runge-Kutta-Merson method) and the absolute 
displacement transmissibility is obtained as 

Td = / ^0 

where |x| is the absolute of the maximum steady-state 

IHclX 

displacement. 

In the absence of viscous damping (C=0) with =0.75, the 
values of so obtained oi^ expressed in dB (-20 log^Q T^) cind are 
indicated in Figure 4.11 by a circle mark along with the results 
obtained by the method of harmonic balance (same as Figure 4.3). 

It can be seen from Figure 4.11 that the numerical results 
also show the stable branch (AC) associated with the anomalous 
jump. To get a clear picture of the response near the anomalous 
jump the phase plots of equation (4.2) are shown in Figures 4.12 
and 4.13 with appropriate initial conditions. Figures 4.12a and 
4.12c show the anomalous response at Q=1.7 and 1,8 respectively 
whereas Figures 4.12b and 4.12d show the resonance responses for 
the same frequencies obtained with a different set of initial 
conditions. Figure 4.13a shows the jump to resonance branch at 
n-2 . 1 obtained with the same set of initial conditions for which 
the anomalous response was obtained in Figures 4.i2a and 4.12c. 
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Figure 4.11 Variation of transmissibility, (in dB) , versus 

frequency, Q. C' = 0 and = 0.75; ( — )stable;( )unstable; 

(o o o) numerical simulation results. 
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Figure 4.14 Phase plot of the ^one-third subharmonic response 
with C=0.0, Cf=0.75, n=4.5 and Aq=Aq= 0 . 01 . 
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Again for Q— 2 , 3 , both the non-resonance and resonance responses 
can be seen in Figures 4.13b and 4.13c. Thus the numerical 
simulations confirm the existence of an anomalous jump in the 
response predicted by the harmonic balance. It was observed that 
the basin of attraction for the anomalous response is very small 
as compared to those of the resonance and non-resonance responses. 

It can be noted from Figure 4.11, that though the response of 
equation (4.2) may not be harmonic, the values of T^, obtained 
from equation (4.28) by numerical simulation match with those 
computed from the assumed harmonic motion. Thus one can rely on 
the method of harmonic balance to provide a good approximation to 
the transmissibility . However, around n=4 . 5 the non-resonance 
response predicted by the harmonic balance differs significantly 
from that obtained by the numerical simulation of equation (4.2). 
A phase plot of equation (4.2) with Q=4 . 5 is shown in Figure 4,14. 
It can be seen from this figure, that this phase plot corresponds 
to a one-third subharmonic in the response which has not been 
accounted for in the harmonic approximation. It is evident that 
the presence of the subharmonic resonance decreases the high 
frequency attenuation rate (given by the slope of the 
non-resonance branch) and consequently the isolator effectiveness. 
The effect of the subharmonic resonances on transmissibility is 
analogous to the wave effects, caused by the isolator inertia, in 
a linear isolator [4]. From the above discussion one can conclude 
that consideration of subharmonic resonances, which occur well 
within the normal operating range of frequencies, is important in 
vibration isolation. 

It is shown in Figure 4.5d that when ^=0 and ^^>71/4, the 
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anomalous jump merges with the main jump. This situation is 
depicted in Figure 4.15 along with the numerical resuits. Again, 
the results obtained by numerical integration match with those 

obtained from 

the harmonic balance method. The phase plots of equation (4.2), 
shown in Figure 4.16, also clearly indicate the merger of the two 
jumps . 

In a soft buffing's oscillator the anomalous response is 

seen to be related with the so-called symmetry-breaking phenomenon 

when the higher order harmonics are included in the solution 

[75-80]. However, in the present case, analysis of the anomalous 

portion including higher harmonics is difficult due to the 

presence of the Coulomb damping and therefore, a numerical 

investigation was carried out. It was seen that in the case of a 

base-excited, hard, buffing's system numerical simulation did not 

reveal any further bifurcations around the anomalous jump region. 

Thus, one may conclude that in the system under consideration, the 
r 

occurence of an anomalous jump in the response is not a signature 
for symmetry-breaking and period-doubling. 

4.3.2 Force excitation t 

The force transmissibility (T^) , given by equation (4.27) 
with calculated from equation (4.23), is shown in Figures 

4.17-4.19 for various values of Cf and C- It may be noted that 
apart from the usual jump observed in a hard buffing equation, no 
anomalous behaviour exists. This can be verified readily following 
the procedure outlined in the Appendix C. In the absence of 
viscous damping (C=0) , the resonance branch extends to infinity as 
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Figure 4.15 Variation of transmissibility, (in dB) , versus 

frequency, f2. C=0 and <^=1.0 ; (-—) stable; ( -) unstable; (o o 

o) numerical simulation results. 


Figure 4.16 Phase plots with <=o and Cf=1.0 (a) n=2.0 and 

Ao=aQ=O.OOl (b) n= 2.0 and Ao=Aq= 3 .o (c) n= 2 .i and Aq=Aq= 0.001 (d) 
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shown in Figures 4.17a-4.17d. However, with a suitable choice of ^ 
one can make the resonance transmissibility bounded as seen in 
Figures 4.18 and 4.19. It is observed from Figures 4.l7a-4.l7d 
that the transmissibility at high frequencies is almost constant 
when the viscous damping is absent. In the presence of viscous 
damping, the high frequency transmissibility is no longer constant 
as can be observed from Figures 4.18 and 4.19. With simultaneous 
increase in both ^ and , the high frequency transmissibility 

increases as shown in Figures 4.17-4.19. 

As mentioned earlier, the approximate method used in the 
present analysis breaks down when . Den Hartog 

[91] presented the exact analysis for a linear restoring force and 
gave the expressions for finite resonance response when C 
Moreover, he also considered the motion with a^ost two stops and 
determined the range of motion without stops. However, the present 
problem with nonlinear spring characteristics is not amenable to 
analytical treatment. 

4.4 Conclusions 

This chapter extends the work of Den Hartog and Ruzicka on 
isolators with combined Coulomb and viscous damping by including a 
nonlinear restoring force. The effect of the nonlinear damping 
terms on the usual jump response of the hard Buffing's equation is 
studied by using the method of harmonic balance. A peculiar 
behaviour in the response of the base excited system is pointed 
out. Analytical results related with the anomalous jump have been 
derived. The existence of an anomalous jump in the 
transmissibility curve of a base-excited, hard, Buffing's 
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isolator, is also confirmed by direct numerical integration. 
Unlike in the case of a soft system, this anomalous jump is not 
associated with the symmetry-breaking and period-doubling. The 
transmissibility indices show that the subharmonic resonances can 
impair the performance of such a nonlinear isolator. The stability 
analysis is carried out and it is found that a transition from a 
stable node to a stable spiral takes place near the break-loose 
frequency. Performance characteristics are presented through a 
detailed parametric study. It is pointed out that with a suitable 
choice of Coulomb damping one can maintain the transmissibility at 
or less than unity even for this nonlinear system. Addition of 
viscous damping contracts the anomalous jump and thereby enhances 
the range of frequency over which the transmissibility can be 
maintained at less than unity. The analysis is done by assuming 
continuous motion and hence it is valid strictly for low damping 
cases. It is known that the secondary resonances and chaos occur 
on the low frequency regime around the jump in the resonance 
response curve of a hard Buffing's equation with linear damping. 
Hence, it is interesting to find out whether or not a suitable 
choice of break-loose frequency (i.e., the Coulomb damping 
coefficient) can eliminate these phenomena. The effects of a 
friction damper on secondary resonances and chaotic motion will be 
taken up in the next chapter. 



CHAPTER 5 


CHAOTIC RESPONSE OF A HARD DUFFING-TYPE VIBRATION ISOLATOR WITH 
COMBINED COULOMB AND VISCOUS DAMPING 


5.1 Introduction 

Numerical simulations of the response of a harmonically 
excited mass on an isolator with a cubic, hard, nonlinear 
restoring force and combined Coulomb and viscous damping are 
presented. These computations have been carried out on an HP-8000 
UNIX system with a subroutine of NAG library using 
Runge-Kutta-Merson method. Appropriate tolerance values have been 
selected to get consistent results. Both force- and base-excited 
systems, in the absence of Coulomb damping, exhibit the 
period-doubling ( with a symmetry-breaking precursor) and (Type I) 
intermittency routes to chaos as the system parameters are varied. 
An addition of a little friction damping quenches the chaotic 
responses in the base-excited system. In the force-excited system, 
however, it is shown that the addition of friction damping does 
not alter the bifurcation structure. 

5.2 Isolation System With Base Excitation 

Let a rigid mass be separated from a harmonically moving base 
through an isolator consisting of a cubic, non-linear spring and 
combined Coulomb and viscous dampers. The equation of motion for 
the mass can be written as (see chapter 4) 
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•* • * ? 0 

A +2CA+2Cf sgn(A) + A = n cos Qx (5.1) 

where the dot denotes differentiation with respect to 
non-dimensional time r. 

5.2.1 Control Of Chaos 

The control of the chaotic responses is recently receiving 
attention of the researchers from several disciplines [107-109]. A 
control strategy for suppressing the chaotic responses in maps is 
presented in reference [108] and has been applied to the Poincare 
map, obtained from an experimental investigation of a 
parametrically driven magnetoelastic ribbon [109]. It is reported 
in reference [109] , that no model for the dynamics is required 

for determining a control strategy. Though several feedback 
strategies have been proposed for vibration control, their 
hardware may be quite expensive and designing a robust controller 
may not be an easy task. Hence, a suitable passive control 
strategy is still an attractive option for a practicing engineer. 
As opposed to the feedback control, the present work examines a 
passive control of chaos based on the conventional vibration 
control methods [3] by incorporating a proper friction damping 
mechanism. 

In fact, the role of the chaotic responses in the design of a 
nonlinear vibration isolation system needs to be critically 
examined. First of all, the amplitudes of the chaotic responses 
are usually quite small compared to those of the resonant 
responses [70]. Secondly, the global bifurcation structure [62, 
68] indicates that the secondary resonances and chaotic responses 
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lie to the left of the primary resonance towards the lo^ 
frequency reqion. Thirdly, the chaotic responses usually occur 
with a high level of excitation and that too within a narrow 
regime of the parameter values. However, in case the chaotic 
responses need to be controlled, it is proposed to achieve this by 
a suitable choice of the break-loose frequency i.e., of the 

friction damping coefficient, To substantiate this suggestion, 

numerical simulation of equation (5.1) is carried out with 
(i.e., only with viscous damping) and C=0-01 by varying Q. As the 
frequency is decreased, the periodic orbit (shown in Figure 5.1a 
for n =0.5) lj!yoses its symmetry due to the symmetry-breaking 
bifurcation and the resulting dual, unsymmetric solutions (at 
Q=0.4) are shown in Figures 5.1b and 5.1c with the corresponding 
initial conditions. Decreasing fi further, a period-doubling 
bifurcation takes place (at Q=0.2) and the resulting period-two 
orbit can be seen in Figure 5. Id . At a still lower value of the 
frequency (n=0.17), a chaotic attractor (not shown here) 
resulting from the usual period-doubling route is obtained. 
Further reduction of shows the emergence of a period-three orbit 
followed by the (Type I) intermittency route to chaos as reported 
in references [53, 70 and 89]. 

Numerical simulation of equation (5.1) for the same value of 
C (=0.01) is now carried out after including a small friction 
damping (C£=0.005). These results are shown in Figures 5.2a to 
5. 2d. Figure 5.2a shows a symmetric, periodic orbit at Q=0.5. It 
can be seen from Figure 5.2b that the symmetry-breaking 
bifurcation occurs again at n=0.4 (same as that in Figure 5.1b). 
Thus, it appears that the symmetry-breaking bifurcation is quite 


Figure 5.1 Phase plots of the period-doubling route with ^=0.01 
and <^=0.0. (a) 0=0.5 and Aq=Aq= 0.1 (b) 0=0.4 and Aq=Aq=0.1 (c) 

0=0.4 and A-=A„= -0.1 (d) 0=0.2 and A =A =0.1. 







-5 


-5 


(c) 

Figure 5.2 


Phase plots with <=0.01 and <^=0.005. 


•(d) 

(a) 0=0,5 and 


Ao=Ao=0.1 (b) n=0.4 and Aq=Aq=0.1 (c) 0=0.2 and Ap=AQ=0.1 (d) 

n=0.15 and Aq=Aj^=0.1 
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insensitive to the value of However, any further reduction of 

Q even up to 0.15 does not reveal the period-doubling cascade as 
confirmed by Figures 5.2c and 5. 2d. Thus, the addition of friction 
damping effectively suppresses the chaotic responses observed in 
the presence of viscous and other non-linear damping forces [89]. 

It may be noted that with the period-doubling occurred 

at around 0=0.2 (Figure 5. Id) and the chaos at Q=0.17. In order to 
eliminate the period-doubling route to chaos, one would tend to 
choose Q^=0.2 (i.e., Cf=(^/8)0^ » 0.0157 [105] ). However, it is 

evident from Figures 5.2a to 5. 2d that, a much lower value of 
(=0.005 with corresponding r2j^=0.113) is sufficient to quench 
chaos. Hence, it can be said that the period-doubling bifurcation 
(unlike the symmetry-breaking) is quite sensitive to the value of 

^f * 

The phase plot of equation (5.1) for $1=0.15, shown in Figure 
5. 2d, indicates the onset of the characteristic stick-slip of a 
friction damper. Such responses can be accurately obtained by 
using numerical integration only after incorporating special 
conditions [94]. However, these studies are not carried out in 
the present work. It has been reported that the low frequency 
chaos and escape phenomenon are common features in soft, Puffing's 
oscillators [44, 47, 52]. The present method of controlling chaos 
with a friction damper is worth exploring in such systems. 

5.3 Isolation System With Force Excitation 

Referring to chapter 4 (see Figure 4.2), if the rigid mass M 
is excited by a harmonic force (f cos wt) and is isolated from an 
immobile foundation, then the equation of motion can be written as 
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K X" + C X' + sgn(x') + Kx^ = f cos wt . (5.2) 

Defining the following non-dimensional parameters . 
r = wt, 

A X = x/Xq , X = x'/{XqU)), X = x"/ix^0)^), < = C/(2Ma)), 

Cj = Cj/(2MXqCJ^) , F = f/(Mw^XQ) with x^ = 


and using them in equation (5.2), one obtains 

X t2CX+2Cf sgn(X) + X^ = F cos z . (5.3) 

In this section, the effect of friction damping on the 
bifurcation structure [62, 68] is examined by analysing the 
results obtained by numerical integration of equation (5.3). 

Equation (5.3) is numerically integrated for various values 
of F with C == 0.05 and = 0.005. The stroboscopic maps (sampled 
at x=2tt) , obtained from these results, are shown in Figures 5.3 
and 5.4. Figures 5.3a to 5.3c indicate the period-doubling route 
to chaos as F is increased from 4.9 to 5.7. Figure 5.3a shows the 
period-doubling , Figure 5.3b shows the period-quadrupling and 
Figure 5.3c represents the chaotic response. Upon further 
increasing the value of F, a period-three solution is obtained at 
F=9.8 (Figure 5.4a). The (Type I) intermittency route to chaos is 
resulted at F=10.0 (Figure 5.4c) with a prechaotic regime (Figure 
5.4b) at F=9.9. When the value of F is increased to 13.4, again 
a period-one solution is obtained (Figure 5.4d). This transition 














(a) 



internuttenoy rout 

to chaos with C-0.05 and Cf=0.005. (a) F=9.8 (b) P=9.9 (c) F=10. 

(d) F«13.4 
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of a period-three orbit to a period-one orbit via a chaotic zone 
has also been observed in references [53, 70]. 

It may be mentioned that Ueda [53] presented a detailed set 
of results for the same parameter values as above and ^^=0.0. A 
comparison of Ueda's results with those presented in this section 
indicates that the bifurcation structure seems to be unaffected by 
friction damping. Numerical simulation carried out for up to a 
value of 0.02 confirms this. It is also observed that the critical 
values of F needed to cause the symmetry-breaking and 
period-doubling are relatively insensitive to the value of 
Just as in the case of a base-excited system, the analysis 
presented in this section does not encompass the stick-slip 
motions . 

5.4 Conclusions 

It has been demonstrated by numerical simulation that the 
instabilities associated with the secondary resonances and chaotic 
motion, encountered in the low frequency regime of a base-excited 
system, can be eliminated by using a suitably designed friction 
damper. The bifurcation structure of a force-excited, hard system 
seems to be unaltered by the addition of friction damping. 


CHAPTER 6 


ROLE OF NONLINEAR DISSIPATION IN SOFT DUFFING’S OSCILLATORS 


6.1 Introduction 

Soft Duffing's oscillators (with a saddle point in the 
unforced situation) are employed as models of various physical ancl 
engineering situations such as Josephson junctions, optical 
bistability, plasma oscillations, buckled beam, ship dynamics, 
vibration isolators and elctrical circuits, etc. , [43-60] . The 
effect of a strictly dissipative force (velocity to the pth-power 
model) on the response and bifurcations of driven, soft Duffing's 
oscillators is considered. The method of harmonic balance is used 
to obtain the steady state harmonic response. An anomalous jump in 
the harmonic response (signifying a break in the resonance curve), 
obtained in the case of linearly damped, soft Duffing's 
oscillators, is shown to persist even in the presence of nonlinear 
damping. It is shown that the bifurcation structure and the 
structure of the chaotic attractor are quite insensitive to the 
damping exponent p. However, the threshold values of the 
parameters, at which bifurcations occur, depend both on the 
damping index and the damping coefficient. The Melinkov criterion 
and an analytical criterion for the period-doubling bifurcation 
have been obtained in the presence of combined linear and cubic 
damping. 
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6.2. Single-Well Potential Oscillator 

Consider the governing equation of motion in the non-dimensional 
form as 


X + 2^ X + 2Cp x|x|^ ^+x+cx^=F cos Qx , c<0,p>0 (6.1) 

where p is the damping exponent, C is the coefficient of viscous 
damping, ^ is the coefficient of pth-power damping, c is the 

Ir 

nonlinearity parameter, F is the amplitude of excitation and Q is 
the frequency of excitation. The dot denotes differentiation with 
respect to time x. In this section we illustrate the effects of 
nonlinear damping on the harmonic response and bifurcation set. 
Assuming the harmonic solution of equation (6.1) in the form 
X = A cos (Qx-cp) (6.2) 

and using the method of harmonic balance one gets the following 
equation for A: 

A^ (1-n^)^ + (9/16) A^ + (3/2) (1-Q^) e a"^ 

+ (2 C A £1 + 2 aP fiP y )^ - = 0 (6.3) 

ir im 

where y is given by 
iP 

yp = (2/v/7T)r[(p+2)/2]/r[(p+3)/2], (6.4) 

r is the standard gamma function. 

It should be noted that equation (6.3) can also be obtained using 
the notion of equivalent viscous damping coefficient [3]. 

Typical curves of A vs. Q with p=2 are shown in Figures 
6. la-6. Id for different values of F. One can see from Figures 
6.1a-6.1c that apart from the usual jump phenomenon, there exists 
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another branch at low frequencies referred to as the 'anomalous 
jump' . When the forcing amplitude F is increased beyond a critical 
value, the anomalous jump merges with the main jump and the 
response curve opens up (similar to the case of an undamped 
system) as shown in Figure 6.1(d). This phenomenon was observed 
with p=l and was related to the symmetry-breaking bifurcation 
[75-80] . 

To obtain the bifurcation set depicting the unstable zone in 
the F-Q plane, one can define 

p = (1-Q^), cr = 2 C and o-p = 2 Cp fiP (6.5) 

and rewrite equation (6.3) as, 

9 + 24 p c + 16 p^ A^ + 16 (o' A + cr A^) ^ -16 F^=0 

Ir 

( 6 . 6 ) 

So long as 'the above equation is bi-cubic in A, for p=l,2 and 

3. In these cases one can employ the condition for repeated roots 
of a cubic equation to get the loci of the turning points (the 
points of vertical tangencies in the response curve) in the F~n 
plane [74]. The bifurcation sets so obtained for p=l and 2 are 
shown in Figures 6.2 and 6.3, respectively. It may be seen that 
the nature of the bifurcation set is similar in both the cases. It 
may be noted that in the soft system considered here (with c<0) , 
there is a limiting amplitude beyond which the spring loses its 
stiffness and the above results are not valid beyond that 
amplitude . 
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6.3 Double-Well Potential Oscillator 

Consider the following non-dimensional equation 

X + 2 X + 2 x|x|P ^ + c(-x+x^) = F cos Qt, c>0, p>0. (6.7) 

In this section we study the effects of on the period-doubling 
bifurcation, for p=l, 2 and 3, using numerical simulation. Also an 
analytical criterion for the period-doubling bifurcation, based on 
the instability of the harmonic solution, is derived for combined 
linear and cubic damping p=3) . Further, the Melinkov 

criterion has also been obtained for such combined damping. 

6.3.1 Numerical simulation 

Equation (6.7) is numerically integrated using Runge-Kutta-Merson 
method to investigate the role of damping exponent on the 
structure of the chaotic attractor. The Poincare maps of the 

chaotic attractors are shown in Figure 6.4 for various values of 
p. It can be seen from Figure 6.4 that the structure of the 

chaotic attractor is insensitive to the value of p. Numerical 

results (not presented here) suggest that both the period-doubling 
and intermittency routes to chaos observed with linear damping 
[38, 89] persist even in the presence of nonlinear damping. 

Though, the bifurcation structure is qualitatively similar for 
various values of p, the threshold values of the parameters at 
which the period-doubling bifurcation occurs depend on the values 
Of both p and Cp. Table 6.1 shows the critical values of F 

required to generate a period 2 orbit with C=0, c=1.0 and n=1.5. 
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Table 6.1 


K? 

1 

2 

3 

0.05 

0.38 

0.38 

0.38 

0.1 

0.41 

0.39 

0.38 

0.15 

0.44 

0.40 

0.39 

0.2 

0 . 49 

0.42 

0.40 

0.25 

0.54 

0.45 

0.41 


It can be observed from the above table that, for the same values 
of C„(>0*05), the critical values of F decrease with increasing p. 
This feature is just opposite to what has been observed in the 
case of a hard, buffing's oscillator (see chapter 3 ) where the 
critical values were seen to increase with increasing p. 

6.3,2 Period-doubling criterion 

A criterion is derived for the period-doubling bifurcation 
based on the instability of harmonic solution predicted by the 
classical stability analysis. Towards this end, we take e=l /2 and 
p=3.0 in equation (6.7). Assuming the harmonic solution of 
equation (6.7) in the form 


X = Aq + A cos {Qx-(l)) 


( 6 . 8 ) 


and using the method of harmonic balance , we get the following 
expressions for A^ and A: 

Aj - (2-3Ah/2 


and (15/8)^ A® - (15/4) (l-n^) A** + 

+ (2CAfJ + 3 C 3 A^bVs)^ - = 


(1-Q^)^ A^ 


0 . 


( 6 . 10 ) 
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To analyse the stability of the solution given by equation (6.8), 
consider Xq=x+ 5 , where 5 is a small perturbation. Employing the 
usual variational procedure [36] one gets the following linear 
differential equation with coefficients varying periodically in 
time: 

5 + cos 2nT] 5 + [Aq+X^cos Qt ^A^ cos 2nr] S = 0 

( 6 . 11 ) 


where Tq = 2C + , 

= -3C3A^n^ 

Aq = 3Aq/2 + 3A^/4 -1/2, (6.12) 

= 3AqA^ 

and A^ = 3A^/4. 

Let us assume 


5= ^ (nT/2)+i?] 


(6.13) 


to get the first approximate solution for 1/2 subharmonic 
instability [78]. Substituting equation (6.13) in equation (6.11) 
and using the method of harmonic balance, one finds that for 
nontrivial solution, the following condition must be satisfied: 


(Aq-qV 4)^ + (C+3C3A2nV2) V -a2/4 = 0 (6.14) 

Equation (6.14) in conjunction with equation (6.10) determines the 
boundary of the period-doubling bifurcation in the F-n plane for 
combined linear and cubic damping. Figure 5 shows the results for 
linear (p=l) and cubic damping (p=3) . These results have been 
obtained as special cases from the analysis presented above by 
assuming, in turn C3=0 and C=0. respectively. It may be pointed 
out that the results for p=l are in perfect agreement with those 
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reported in reference [78]. It can be seen from Figure 6.5 that 
the excitation amplitude needed to generate the period-doubling 
bifurcation for p=3 is lower than that for p=l. So this trend is 
in conformity with the numerical results presented in Table 6.1. 

6.3.3 Melinkov Criterion 

The Melinkov criterion for the double-well potential 
oscillator has been obtained in the case of linear damping [43]. 
It has been pointed out that the Melinkov criterion depicts the 
occurence of fractal basin boundaries rather than the onset of 
chaos [44]. Following reference [43], the Melinkov criterion for 
combined linear and cubic damping has been obtained for c=0.5 as 

F = (| C + 8 C 3 /"sec'^ht tan'^ht dt) (6.15) 

““00 

Evaluating the integral in equation (6.15) one gets, 

^ = <1 c - ii = 3 ) 

The coefficient of C 3 iri the above equation is smaller than that 
of C- One can, thus, conclude that the fractal basin boundaries, 
just like the period-doubling bifurcation, occur for a lower value 
of F in the case of cubic damping as compared to linear damping 
(of course with ^=^ 3 ) • 

It should be noted that the present analysis can be easily 
extended to the case of a harmonically driven pendulum (with 
periodic potential) [76]. Following reference [50] it can be shown 
that for a driven pendulum the Melinkov criterion is obtained as 
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F 


(8 c + -r C3) 


cos h(nQ/2} 
n 


(6.17) 


It is worthwhile to note the difference in the coefficients of ( 
and ^2 equations (6.16) and (6.17) . For the driven pendulum 
(which is also a soft Puffing's oscillator) , the critical value of 
F, predicted by the Melinkov criterion, is more for p=3 than that 
for p=l with C=C 3 ■ This is in contradiction to what has been 
just observed for the double-well potential oscillator. 


6.4 Conclusions 

A detailed investigation of the effects of strictly 
dissipative nonlinear damping on the harmonic response, 
bifurcation set and chaotic motion of harmonically driven soft 
Puffing's oscillators is carried out. Numerical results indicate 
that the structure of the chaotic attractor and the routes to 
chaos are quite insensitive to the value of the damping exponent 
p. A parametric study is reported to show the influence of the 
nonlinear damping on the onset of period-doubling route to chaos. 
An analytical criterion has been derived, in the presence of 
combined linear and cubic damping, to estimate the critical values 
of the excitation at which the period-doubling bifurcation takes 
place. The Melinkov criterion for such systems has also been 
obtained. These results indicate that the correct model of damping 
nonlinearity is important in accurate prediction of the onset of 


chaos. 


CHAPTER 7 


EXPERIMENTAL INVESTIGATIONS 


7. 1 Introduction 

In this chapter, some preliminary experimental results 
are presented for a nonlinear vibration isolator made in the form 
of a hollow rubber tube. This isolator has softening type 
restoring force characteristics. The jump phenomenon has been 
observed in the experiments. The transmissibility characteristics 
of this system obtained experimentally indicate that the soft 
characteristics of the hollow rubber tube are beneficial for 
effective vibration isolation. 

7.2 Experiment 

The isolation system under consideration is shown in 
Figure 7.1. A hollow rubber tube (length=102mm, inner 
diameter=24mm, outer diameter=32mm and Hardness=40 IRHD) clamped 
to the projections in the top and bottom plates acts as the 
isolator. In addition to the internal damping of the rubber, 
provision is made to include two rubber pads with adjustable 
pressure which can provide additional damping as the top plate 
rubs against them. The static stiffness of the isolator is 
obtained from a compressive test carried out on a universal 
testing machine. The resulting force-displacement curve, shown in 
Figure 7.2, indicates a soft characteristic of the isolator at 
small deformation. However, as the load is increased further, the 
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Damping pads 
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Figure 7.1 Experimental set-up 
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Figure 7 . 2 


Stiffness characteristics 
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isolator has an almost linear stiffness characteristic. 

In order to measure the transroissibility characteristics, the 
bottom plate of the isolator is clamped on to a mechanical shaker 
which provides a sinusoidal base excitation whose frequency and 
amplitude can be set to desired values. A mass (1,95 Kg) is fixed 
to the top plate of the isolator. The vibration of the mass is 
sensed by an accelerometer (B&K 4370), and' the output signal is 
fed to a Digital storage oscilloscope and a H.P. real time 
spectrum analyser via a charge amplifier (B&K 2635) (see Figure 
7.1). 

The displacement transroissibility (T^) is obtained as the 
ratio of the maximum displacement of the mass to the amplitude of 
the base excitation. The transroissibility values are obtained for 
both increasing and decreasing values of frequency. The results 
obtained are shown in Figures 7.3 and 7.4 for two different values 
of the base amplitude. The transroissibility curves are bent to the 
left clearly indicating the softening nature of the isolator. To 
depict the jump phenomenon clearly, the zone around the resonance 
is expanded as shown in Figures 7.5 and 7.6. It may be observed 
that as the base excitation is increased the jump region also 
increases as expected. Also the value of the frequency at which 
the maximum transmissibility occurs is shifted to a lower value 
(from 33.2 Hz to 28 Hz) and the value of the maximum 
transmissibility increases as the base excitation is increased. 

A curious observation worth mentioning is that the transients 
are very long at the frequency when the jump occurs. This 
phenomenon is analogous to the case of a first order phase 
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Figure 7.3 Variation of transmissibility with frequency with base 
amplitude = 0.28 mm. x x x Increasing frequency; o o o Decreasing 
frequency ; 



(QP) Pi 


Frequency (Hz) 


Figure 7.4 Variation of transmissibility with frequency with base 
amplitude = 0.625 mm. x x x Increasing frequency; o o o Decreasing 
frequency; 



Td (dB) 


4 

Frequency (Hz) 

Figure 7.5 Jump phenomena with base amplitude = 0.28 mm. x x 
Increasing frequency; o o o Decreasing frequency; 



Figure 7.6 Jump phenomena with base amplitude - 0.625 mm 
Increasing frequency; o o o Decreasing frequency; 


XXX 
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transition. It appears that the role of noise on the onset of 
bifurcation is the cause for these long transients. Another point 
to note is that the values of obtained while decreasing the 
frequency are found to be more than those obtained while 
increasing frequency around the jump region. This could be 
attributed to the fact that due to the long transients mentioned 
above, there is a considerable increase in the temparature of the 
rubber isolator (with associated change in the property values) 
while decreasing frequency. 

One important observation from the view point of vibration 

isolation is, that as the base amplitude is increased (Figures 7.3 

and 7.4), the frequency at which (in dB) becomes negative 

decreases substantially from 60 Hz to 49 Hz, This indicates that 

there is a substantial increase in both the operational frequency 

range and the high frequency attennuation rate due to the 

softening nature of the isolator. However, as mentioned above the 

maximum value of T, is increased from 14.8 dB to 16.2 dB as the 

d 

base excitation is increased. This increase in T^ can be offset by 
suitably increasing the damping. To see the effect of increasing 
the damping, two settings of the adjustable damping pads are used 
and the resulting transmissibility values are plotted in Figures 
7.7 and 7.8. For both these settings of damping, the base 
amplitude is fixed at 0.625 mm. Comparing Figures 7.4, 7.7 and 7.8 
it can be said that the increased damping eliminates the jump, 
reduces the maximum value of T^ and increases the frequency at 
which T^ (in dB) becomes negative. As expected, the high frequency 
attennuation rate is poorer with higher damping. But a suitable 


(ap) Pi (gp) PjL 
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Figure 7.7 Variation of transmissibility with frequency with base 
amplitude = 0.625mm and increase in damping, x x x Increasing 
frequency; o o o Decreasing frequency; 
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Frequency (Hz) 


Figure 7.8 Variation of transmissibility with frequency with base 
amplitude = 0.625 mm and further increase in damping, x x x 
Increasing frequency; o o o Decreasing frequency; 
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combination of a soft nonlinear spring and the optimum value of 
damping can yield a better design than the linear isolator. 

The Fourier spectrum of the response of the isolator is 
obtained by feeding the signal to a real time spectrum analyser. 
As the base amplitude is increased, a 3/4 subharmonic appears in 
the response as shown in Figure 7-9. At some other excitation 

I 

frequency with the same base amplitude, a 1/2 subharmonic is also 
revealed as shown in Figure 7-10. However, no further bifurcations 
could be observed with this level of excitation. 





CHAPTER 8 


CONCLUSIONS 


In the present thesis, harmonic and chaotic responses of a 

I 

single degree-of-freedom vibration isolation system with 
nonlinearity in both stiffness and damping have been studied. 
While the method of harmonic balance has been used to obtain the 
harmonic response, the chaotic responses are investigated through 
numerical simulation. Both force- and base-excited systems have 
been considered. Also special attention is paid to understand the 
role of strictly dissipative nonlinear damping in the chaotic 
dynamics of such systems. The main conclusions that can be drawn 
from this work, are summarised below: 

1. Detailed performance characteristics of nonlinear vibration 
isolators presented as transmissibility vs frequency plots suggest 
that, an isolator with a soft characteristic is superior to that 
with a hard characteristic. With increasing damping index p, the 
jump width reduces and consequently the unstable zone becomes 
narrower. It is seen that, for the parameter values considered, an 
isolator with an asymmetric restoring force may not perform 
satisfactorily for base excitation. 

2. It is shown that both the period-doubling and intermittency 
routes to chaos occur for all values of p, thus preserving the 
bifurcation structure of a linearly damped, Duffing's oscillator. 
It is pointed out that nonlinear damping can be effectively used 
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for controlling chaos, 

3. An anomalous jump in the response of a base-excited/ hard 
Duffing^type vibartion isolator with combined Coulomb and viscous 
damping is shown to exist when the method of harmonic balance is 
used. This is also confirmed by direct numerical integration. 

4. It is shown that with a suitable choice of Coulomb damping one 
can maintain the transraissibility at or less than unity over the 
entire frequency range even for a nonlinear system under 
base-excitation . 

5. It has been demonstrated by numerical simulation, that the 
instabilities, associated with the secondary resonances and 
chaotic motion, encountered in the low frequency regime can be 
eliminated by using a suitably designed friction damper. The 
bifurcation structure of a force-excited, hard system seems to be 
unaltered by the addition of friction damping. 

6. With the same value of the damping coefficient, the threshold 
values of forcing required to cause period-doubling bifurcation 
are seen to decrease with increasing p in the case of a 
double-well potential oscillator. This feature is just opposite to 
what has been observed in the case of a hard, Puffing's oscillator 
where the critical excitation values are seen to increase with 
increasing p. 

7. Experimental results indicate that a suitable combination of a 
soft nonlJLnear spring and the optimum value of the damping can 
yield a better design than a linear isolator. 
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8.1 Suggestions for future work 

The following directions may be taken up to extend the work 
carried out in the present thesis: 

1. The stability analysis of isolators with Coulomb damping should 
be expanded to include the stick-slip phenomena. Also more 
realistic models of Coulomb damping with variable friction force 
need to be considered. 

2. The emphasis of the experiment carried out in the present work 
is to obtain a qualitative understanding of a soft isolator 
subjected to a base excitation. No attempt is made to model the 
isolator and obtain the quant itaUve predictions. A detailed 
modelling of the damping phenomena needs to be carried out. Also 
the design of friction damper to exploit the break-loose phenomena 
in the case of base excitation needs to be examined. 

3. Characterisation of the chaotic attractor via Liapunov 
exponents, fractal dimensions and entropies also deserves 
attention. 
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APPENDIX A 


STABILITY ANALYSIS OF A NONLINEARLY DAMPED DUFFING’S OSCILLATOR 


A. 1 Introduction 

The equation of motion of a harmonically excited, hard, 
Duffing 's oscillator with generalized pth-power damping model 
[ 112 ] , is given by 

X + C' X|X|P~^ + aX + I3X^ = F' cos u' t (A.l) 


where a>0, /3>0, p>0 and the dot denotes differentiation with 
respect to time t. The parameters C and p are called the damping 
coefficient and exponent, respectively. The harmonic solution of 
equation (A.l) is sought in the form 

X^ = A' cos (w't “ (p) . (A. 2) 
which can be easily obtained by any standard method like the 
equivalent linearisation technique or the harmonic balance method. 

In the standard variational technique [36] of stability 
analysis of the solution given by equation (A. 2), a small 
variation t) is considered in the form 


Substituting equation (A. 3) in equation (A.l) and using the fact 

that X satisfies equation (A.l), one gets a linear differential 
o 

equation in t) with time varying coefficients (after neglecting 
terms higher than linear in tj) . This approach has been termed 
dynamic analysis [113] as opposed to the static analysis where the 
bifurcation of the steady state harmonic solution is investigated 
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through the turning points heving vertical tangents to the 

response curves in the A'-o)' plane. Only the dynamic analysis can 

capture unstable regions of secondary and higher orders. The 

dynamic analysis is increasingly becoming popular towards 

obtaining approximate analytical criteria for the onset of 

symmetry-breaking and chaotic motion of a Duffing 's oscillator 

% 

[114]. 

It must be emphasized that the dynamic analysis cannot 
be carried out in the usual manner unless the damping exponent p 
is an odd integer. Moreover, in the range 0 < p < 1, the 

derivative of the damping force is discontinuous. Consequently, 
approximate methods have been attempted in the past. For example, 
in reference [51], quadratic damping (p=2) has been replaced by a 
linear plus cubic model. Some other methods have been proposed 
while carrying out the stability analysis of the responses of 
impact dampers [115] and hysteretic oscillators [ill]. 

In this appendix, another approximate method is proposed 
using the notion of equivalent viscous damping coefficient. The 
method is based on the correspondence between the static analysis 
and the primary unstable region obtained from the dynamic analysis 
for p - 1 and 3. Following the usual method of variational 
technique, one obtains a Mathieu equation for p = l and a Hill 
equation for p = 3 . It is shown that one can get away with only 
Mathieu equation for all values of p by suitably modifying the 
equivalent viscous damping coefficient in this variational 
equation. This simplification has been justified from the fact 
that strictly dissipative, velocity-dependent nonlinear damping 
alone does not give rise to any instability. Numerical results 
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for p = 2 obtained by the present method show good agreement with 
those obtained by another approximate method [111] which uses, 
after some approximation, the van der Pol plane. Unlike this 
van der Pol plane method, the present analysis, however, can also 
be used to obtain approximately the higher order unstable regions. 


A. 2. Static Analysis 

The static analysis or the classical jump phenomena in the 
response curve for the harmonic solution of equation (A.l) is 
carried out as follows. It is known that for the pth-power 
damping model, the equivalent viscous damping coefficient can be 
obtained as [112] 


C' 

e 


C' A' 


to' y 


P 


(A. 4) 


where 


^ _ _2 rr(p+ 2 )/ 2 ] 

^r[(p+3)/2] 


(A. 5) 


with r representing the standard gamma function, 
Now equation (A.l) can be rewritten as 


X + C' X + aX + IBX = F'cos Lo't 

% 


(A. 6) 


and the amplitude A' of equation (A. 2) is known to satisfy the 
relation [3] 


A 


,2 _ 




C'^ + 

% 


to'^ - a 


>'2 


2 


(A. 7) 
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Introducing the following non-dimensional parameters: 

(2p-3)/2 -j 1/t 

OJ = coVv/a, Cp = C' JI^ZTY72' ^ = (/3/«) A', F={IB/a) 

and C = C (A. 8) 

e„ P P 


one can write equation (A. 7) as 


A 


2 




1 - 


3 a2>2 


(A. 9) 


It has been shown [35] that irrespective of the value of p, the 
response curve given by equation (A. 9) yields the usual jump 
phenomenon observed in the case of p=l. The use of equivalent 
viscous damping coefficient is nothing but the method of 
equivalent linearisation [3]. The same result, as given by 
equation (A. 9), is also obtained by the method of harmonic balance 
[35]. 

For a given value of F' , differentiating both sides of 
equations (A. 7) with respect to w' , one obtains an expression for 

dA ^ ♦ I 

• Setting the denominator of this expression to zero, i.e, 

dA' 

^ CO, one gets the locus of the points in the A'-w' plane 
having vertical tangents to the response curves. The expression 
thus resulted can be put in the following non-dimensional form; 


p w^+ 
"P 



1-3 

4 





1-3 A^ 
4 


0 


(A. 10) 
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One can rewrite equation (A. 10) as 


1 

2 

r 2-, 

3 A^ 

2 


L 0,2 J 


4 2 

L (j J 



(A. 11) 


Introducing a new set of non-dimensional parameters 5 and e 

given by 


Ce 1 9 

2^ ^ a ^ and c = I ^ 


(A, 12) 


equation (A. 10) can be recast in the form 


9 9 9 

(5 - 1 )^ = 

or, 5 = 1 +/ . (A. 13) 

ir 

Attention may be drawn to the appearance of the term v^p in the 

parameter fi and dependence of C on A, u) and p. 

^ ®P 

A. 3 Dynamic Analysis 
A. 3.1 Linear damping (p=l) 

With p = 1.0, i.e., in the case of a viscously damped 

Duffing' s oscillator, the variational technique yields the 
standard damped Mathieu equation in rj [36], 

2 

^ + [5 + 2 c cos 2z] = 0 

dz 


(A. 14) 
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where jLi^, 5 and c are given eguation (A. 12) after substituting p = 

1 (i.e., C = C, = ) and z = w't - 0. 

' e 1 

P 

The transition curves (in the 5-e plane) bounding the primary 
unstable zone of equation (A. 14) are given approximately by 
[36-40] 


5 = 1 ± 


c - 4M 


2l 2 
^ ' + 0(e"^) 


(A. 15) 


which is nothing but equation (A. 13), thus confirming the fact 
that the loci of the points vertical tangencies of the response 
curve in the amplitude-frequency plane correspond to the 
boundaries of the primary unstable region of equation (A. 14). It 
is easy to show that the curves bounding the primary unstable zone 
are given on the A-co plane by [112] 



a2 8 , 2 , , 4 

A = _ (oj -1) + - 

[(w^-1)^ - 

2 2 ^^'^ 

3C^w^] 

(A. 16) 

and 

a2 8,2,, 4 

A = g (u -1) - - 

[(m^-l)^ - 

2 P 1/2 

3C^W^] 

(A. 17) 

The 

transition curves 

bounding 

the secondary unstable 

zone of 


equation (A. 14) are given by [36-40], 

1/2 

5 = 4 + I + (cVi6-16m/) + 0 (c^) . (A. 18) 


In the absence of damping (i.e., with p = 0 in equation 
(A. 18)), employing equation (A. 12) the transition curves bounding 
the secondary unstable zones can be obtained in the A— w plane as 
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I 



and 




3 - 11/2 
2 


A 


2 



15 

160 )^ 



(A. 19) 

(A. 20) 


Following the results given in reference [37] , it may be observed 
that eguation (A. 19) is guite accurate for w > 0.5 whereas 
equation (A. 20) is accurate only in the range 0.5 < cj < 0.8. 
Figure A.l shows the primary and secondary unstable zones in the 
A-cj plane for the undamped Duffing's oscillator. In this figure, 
curves I (the backbone curve) and II (the locus of the single 
turning points) are given by equations (A. 16) and (A. 17) 
respectively with = 0. It can easily be shown from equations 
(A. 16), (A. 17) and (A. 12) that for this undamped situation, the 
primary unstable region bounded by curves I & II is quite accurate 
for the entire range of frequency. The secondary unstable region 
starts from w = 0.5 and equation (A. 19) produces the curve CPQ 
which, beyond the point Q, deviates drastically from the exact 
curve CPQ' . The exact curve can be obtained through numerical 
computation of the transition curve of the Mathieu equation. In 
the same figure, the response curve for F = 1.0 is also shown 
whose unstable portiphs are indicated by the dashed lines. 

The secondary unstable region in the presence of viscous 
damping can be obtained through the use of equation (A. 18). The 
two transition curves can be put together as 





Figure A.l Primary and Secondary instability regions with 
typical response curve for an undamped Buffing's oscillator 






161 


Now, substituting for S, c and jj.^ from equation (A. 12) (with p = 

1) in the above equation and simplifying one gets the following 
. . 2 

quartic equation in A : 

a^(A^)^ + a^CA^)^ + a^CA^)^ + a^(A^) + a^ = 0 (A. 22) 

where a^ = 45/256 , a^ = | , a^ = - 48w^, a^ = 192 (j® - 48w'^, 

and ap. = 128a)® - 160)"^ - 256 o)® - 64C^ o)®. 

5 e3_ 

Solving equation (A. 22) numerically and retaining only the real, 
positive roots one can obtain the secondary unstable region in the 
A-o) plane. Figures A. 2 and A. 3 show respectively, the primary and 
secondary unstable regions for some low values of damping. The 
representation of the primary and secondary unstable zones on the 
A-o) plane may help in locating the chaotic response of equation 
(A.l) with reference to the classical jump curves [114]. 


A. 3. 2 Cubic damping (p=3) 

With p=3 in equation (A.l), considering the variation of 
the solution given by equation (A. 2) and carrying out the usual 
procedure, one gets the following non-dimensional equation: 


j2 r4ji_ 

dz L /3 


4]U 

cos 2z 

y/3 


d7) 

dz 


+ [5 + 2c cos 2z] T 7 = 0 


(A. 23) 


Where S and e are still given by equation (A. 12) with p=3 and 

3 

Z = u't-(p. 

The Lindstedt - Poincare technique [38] can be used to 
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Figure A. 3 Secondary unstable regions for various coefficients of 
linear damping. 
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obtain the tra-nsition curves for the primary unstable zone. 

A 

Towards this end, one sets = c and writes 


7](z, c) = v^iz) + c 17 .( 2 ) + e 17 ,( 2 ) + 


(A. 24) 


and 5(c) = 5^ + c8^ + e 5^ + 


(A. 25) 


Substituting equations (A. 24) and (A. 25) in equation (A. 23) and 
equating coefficients of like powers of e, one obtains 


d 17 , 


dz 


2^ + ^o ^o = ° 


(A. 26) 


d T). 


dz 




■4U. 


A 

4h, 


1 dT), 


L /3 ]/3 


cos 2z 


dz 


[6^ + 2cos 2z] T7 q (A. 27) 


and 


^2 

d 17 , 


dz 


r- + «o ^2 


4jU, 


1 di). 


L /3 


COS 2 z 


dz 


■^2 ”0 


[5^ + 2cos 2z] Tj^ 


(A. 28) 


According to the Floquet theory [36], -q has the period either tt or 
271 along the transition curves. 

Therefore, one can write 


■>7q = cos nz + sin nz 


(A. 29) 


where n is a non-zero integer related to So by 


n = 5 




(A. 30) 



and and are arbitrary constants at this point. 

For the transition curves bounding the primary unstable zone, 
unity and hence 

= 1 and = a^ cos 2 + b^ sin z. 
substituting for into equation (A. 27) and simplifying 

obtains 


d^T?, 


dz 


r" + "o ^1 = 


A 

■4/i. 


L /3 


i - b 5, + i f 
o o 1 2 ( 


A 

4W- 

— — a^ + 2b„ 
>/3 ° ° 


Sin z 


4^3 ^ 

— — ' b -f a. 5 - -ir 
L v/3 o 0 12 


/4|U- 


— ^ b “ 2a 
/3 o 


COS z 


A 

1 ) 1 (^^3 

' a_ + 2b sin 3z + i — — b - 2a_ | cos 3z. 

oj 2(^0 


✓3 ° 




Thus the 
given by 


conditions for eliminating secular terms 


A 



211 ^ N 


a + 
o 




0 


from T)^ 

(A 


and 


(1 + 5^) a^ + 





(A 


For a nontrivial solution to exist for a^ and b^, one equates 
coefficient determinant for equations (A. 32) and (A. 33) to 
and obtains the following equation : 


— + (1"'4 IJ.2 ) 


165 

n is 

one 


.31) 

are 

.32) 

.33) 

the 

zero 


1/2 


(A. 34) 
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Substituting for 5^ and 5^ in equation (A. 25), the transition 
curves for the primary unstable zone are obtained as 

1/2 

5 = 1 + (c^ - 4 /i^) +0 (c^) (A. 35) 


Once again it can be seen that equation (A. 35) is nothing but 
equation (A. 13). 

For deriving the transition curves bounding the secondary unstable 
zone it is assumed that iJ.^= O(e^) = carrying out a 

similar procedure outlined above, one gets 


5 




(A. 36) 


A. 3. 3 Proposed method 

Comparing equations (A. 15) and (A. 35) one can conclude 
that the form of the expression defining the bounds of the primary 
unstable zone is identical for p = 1 and p =3. The only difference 
is in the corresponding subscript used for /i. It should be noted 
that equation (A. 23) differs from equation (A. 14) since the 
coefficient of ^ in the former varies periodically with time. 
However, as shown in the appendix B, the time variations of the 
coefficient of alone, in equation (A. 23), does not result in 
any instability. Therefore, with time varying coefficient of v, 
it is proposed to approximate the time varying coefficient of ^ 
in equation (A. 23) by some average value. Comparing equations 
(A. 15) and (A. 35) it is clearly seen that for p = 3, this average 
value turns out to be V3/2 times the constant part of the 
coefficient. This approximation may be tested by obtaining 
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results for the secondary unstable zones. The transition curves 
bounding the secondary unstable zone are given by equation (A. 36). 
Equation (A. 18) with replaced by yields the approximation of 
this secondary unstable zone. a comparison of equation (A. 36) 
with equation (A. 18) indicates that the approximation, though a 
little more conservative (resulting in a wider secondary unstable 
zone ) ! is quite good, specially for low values of the damping 
coefficient. 

As mentioned in section i, for a general value of the 
damping exponent p, the variational method is not amenable for 
analytical treatment. Keeping in mind that equations (A. 15) and 
(A. 18) give satisfactory results even for p=3 when is replaced 
by as discussed in the proceeding paragraph, it is proposed to 
use the same equations for any value of p (o<p<3) with replaced 
by Mp. The proposed method is illustrated for the following set 
of parameter values ; p=2.0, F'=1.0, o)'-2.0 and C'=0.O5. The 

three roots of the amplitude of the response can be obtained by 
using equation (A. 8) as 

A^ = 0.34339, A 2 = 1.847, A^ = 2.1022. 

The first root corresponds to the non-resonant branch (stable) , 
the second to the unstable branch and the third one to the 
(stable) resonant branch. The stability/instability of these three 
roots are confirmed by the van der Pol plane method outlined in 
reference [ill]. In this reference also, an approximation was 
made by neglecting the higher harmonics in the Fourier expansion 
of the damping force. Using equation (A. 12) the value of S, c and 
corresponding to the above three roots of A are obtained as 

A 

given below: 
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0.0103 
0.0554 
A^ 0.0631 


5 

0.2942 

1.5293 

1.9073 


c 

0.0221 

0.6397 

0.8286 


The transition curves bounding the primary unstable 

zones for p=2 given by equation (A. 15) with /i^ replaced by 
shown in Figures A,4a-A.4c. The shaded portions correspond to the 
unstable regions. The points corresponding to the roots of A are 
also indicated in the some figure by a cross mark. One can see 
that points A^ and A^ lie outside the shaded portions as shown in 
Figures A. 4a and A. 4c, respectively, and hence stable; where as 

the point A 2 is inside the shaded portion as shown in Figure A. 4b 

and hence unstable. Thus, the results are consistent with those 
obtained by using the method of reference [111]. Similar 

numerical checks have confirmed the validity of the present 
approximate method of obtaining the primary unstable zone by using 
equation (A. 15) for a number of values of p in the range 0 < p 
3. The method, however, cannot be used for p=0 i.e., a Duffing's 
oscillator with Coulomb damping. 


A. 4 Conclusions 

An approximate procedure using a modified equivalent 
viscous damping coefficient and the variational formulation, has 
been proposed to obtain the primary (and higher order) unstable 
zones of the steady state harmonic solutions of the Duffing's 
oscillator with pth-power damping model. Simple Mathieu equation 
is used irrespective of the damping exponent p. Numerical results 
obtained for the primary unstable zone are found to be in 
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Figure A. 4 Primary unstable zones for quadratic damping (p=2) 
^2 = 0.0103; (b) jU 2 = 0.0554; (c) n^=0.Q631 



agreement with that predicted by another approximate method. An 
analytical procedure to depict the primary and secondary unstable 
zones on the amplitude-frequency plane is outlined which may 
enable one to obtain approximate analytical criteria for chaotic 
responses based on the stability of the harmonic solution. 


APPENDIX B 


STABILITY ANALYSIS OF AN OSCILLATOR WITH LINEAR STIFFNESS AND 

CUBIC DAMPING 


Consider an equation of the following form [116] 

ii + d(t) 71 + e(t) 7 ) = 0 (B.l) 

where d(t) and e(t) are continuous and differentiable functions of 
time. Following reference [116], one obtains the conditions for 
stability of equation (B.l) as 

e(t)>0 (B*2) 

and 2 d(t) e(t) + e(t) £ 0 . (B.3) 


Now consider a system with a linear spring and a cubic 
damper (by putting p=3 and |3=0 in equation (A.l)) governed by 

X + C'X^ + a X = F' cos w't with C' > 0 and a > 0 (B.4) 


Applying the variational technique to the harmonic 
solution of equation (B.4), one obtains. 


4* 

71 + 


^^3 

/3 


4^3 

COS 2Z 


71+ 71-0 


with z=a)'t-0 


(B.5) 


Comparing equation (B.5) and (B.l) one can write that 


171 


172 


d(t) = i- [ 1-cos 2z] 

\/3 

and e(t) = 1.0 . (B.6) 

Noting that in equations (B.6), d(t) ^ 0, e(t) > 0 and e(t) = o, 
one can easily see that the conditions for stability of equation 
(B.5), as dictated by equations (B.2) and (B.3), are satisfied. 
This implies that there are no unstable zohes for equation (B.5) 
and consequently there is no jump phenomenon in the response of 
equation (B.4) [3]. This analysis shows that the time variation 

in the coefficient of Tr| alone in equation (B.5) does not give rise 
to any unstable zone. 



APPENDIX C 


SOME ANALYTICAL CRITERIA 

In this appendix, some simple analytical expressions for 

quantities related to the anomalous jump response have been 

derived. The response equation, for the base excited system given 

by equation (4.6), can be put in the form of a cubic equation with 

2 

the transformation , to yield: 

+ a^ A^^ + 32 + 33 = 0 (C.l) 

where 3 ^^= -2iQ^/ a, and 33 = (4 /3^ - Q‘^)/a^ . (C.2) 

Equation (C.l) can be put in the standard form with the 
transformation A 2 = A 3 - a^/3, resulting in the following equation 

A3^ + q^ A3 + g2 = 0 (C- 3 ) 

2 

Where q^^ = - a^ /3 + 

(C.4) 

and q 2 = 2 27 - a^ 33 / 3 + 33 . 

Applying the condition for repeated roots of equation (C.3) [120], 
we get 

q 2 ^/ 4 + 27 = 0 ==> 4/3^ = 0 (C.5a) 

or 4^^ Cf^ - t 4/ (27a) = 0 . (C.5b) 
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The first of the above equations corresponds to the 

. 1/2 

break-loose frequency. For Q > (i.e., Q > (S^^/rr) ), the 

values of n obtained by solving equation (C.5b) for a particular 

correspond to the critical frequency of the anomalous jump (Q ) 

and the non-resonance critical frequency of the usual jump (n^) . 

The values of Q and n , so obtained, are shown in Figure 4.4 for 

an 

various values of 

Now it is interesting to find out the value of = ^^f^m 

which n and Q coincide i.e., when the anomalous jump just meets 

the usual jump. This corresponds to the condition for repeated 

roots of equation (C.5b). Noticing that equation (C.5b) is also a 
• . . 2 

cubic equation in Q ( say W) and substituting for a and (3 in 
(C.5b), one obtains 

3 442 42 22 

- 3"*/ 2^^ + 3^ 2V(7T)^ Cf == 0. (C.6) 

The condition for repeated roots of equation (C.6), is obtained, 
in case > 0, as 

(2Cf/’T)^ = 3^/2^° => (Cf)j^= TT w 71/4. (C.7) 

Beyond this value of ( i.e., > (Cf)j^)/ the usual jump opens 

up as in the soft Buffing's case [77], and merges with the 
anomalous jump. The break-loose frequency corresponding to (Cf)j^ 
is given by (see equation (4.8)) 


(C.8) 




